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SUMMARY 


The first part of the present report deals with the theory of 
elasticity of orthotropic plates. The general differential equation 
of the bending surface of orthotropic rectangular plates whose principal 
directions of stiffness are not parallel to the plate edges (general— 
orthotropic plate) is indicated. The possibility of applying the theory 
of elasticity of orthotropic plates to plywood is investigated and the 
relationship between the elastic moduli and the angle between the load 
direction and the principal stiffness directions is described for 
plywood of various constructions. 

The second part is concerned with the stability equation of the 
general— orthotropic rectangular plate under uniform shear and axial 
loads. This differential equation, which, compared to the differential 
equation of the orthotropic rectangular plate whose principal directions 
of stiffness are parallel to the plate edges, contains additional terms, 
can be solved for a very long plate strip by the Southwell— Shan formula 
for the isotropic plate strip. Besides the exact solution an approximate 
solution is given. For the most important practical case of orthotropic 
plate strips, whose principal directions of stiffness are inclined at k^>° 
with respect to the plate edges, the buckling loads and buckling lengths 
in pure compression and pure shear are plotted against the stiffness 
values of the plate. For the special stiffness values of plywood plates 
of various constructions, the buckling loads and buckling wave lengths 
were determined for all inclinations of principal stiffness directions 
( 0 ° to 180 °) relative to plate edges and also plotted for pure compressive 
and pure shear loading. 


*"Beitrag zur Frage der Beulung orthotroper Platten, insbesondere 
von Sperrholzplatten . " 
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I. INTRODUCTION 


Thin plyvood strips are frequently used in airplane designs as spar 
wets, wing covering, etc. They are subjected to forces within the total 
assembly of the structural member which usually lie in its plane and 
stress the plates in compression or shear. Since the plates are fre- 
quently cut so thin that buckling under such loads is to be expected, 
the knowledge of the critical loads at which the plates change from 
their originally straight, stable position into a new, buckled position 
should be of interest. 

Being built up of plies at right angles to each other, the plywood 
plates are regarded as orthogonally-anisotropic plates (orthotropic 
plates according to the term introduced by Huber). By its construction 
the stiffness of the plate is dependent on the position of the elastic 
axis relative to the direction of the fiber (grain) . 

Although the plywood plate cannot be termed completely homogeneous 
because of the violent property flunctuations over the wall thickness, 
the subsequent investigations of the buckling stiffness of orthotropic 
plates nevertheless are generally made for a homogeneous orthotropic 
and elastic material. The effect of the incomplete homogeneity of the 
plywood plate on the stiffness, etc., is particularly investigated. 

The problem of stability of orthotropic plates has been treated 
by C. Schmieden (reference 1), Bergmann-Reissner (reference 2), and 
Seydel (reference 3), with special reference to the case of pure shear 
loading in connection with the Southwell-Skan investigations on the 
stability of isotropic plates (reference 4). Basic investigations on 
the orthotropic plate have been made by Huber (reference 5 ) who also 
set up the differential equation of bending of the orthotropic plate. 
Further data on the stability of orthotropic plates, including the 
data on buckling load under compressive stress, are given by Timoshenko 
(reference 6). 

But all these reports deal only with orthotropic plates with 
principal directions of stiffness (on plywood, the directions of the 
grain of the inner and outer ply) parallel to the edges of the rectan- 
gular plate (termed "special orthotropic plate" hereinafter). For the 
use of orthotropic plates, especially of plywood, as spar or wing 
covering, however, the knowledge of the buckling strength of plates 
with principal stiffness other than parallel to the edges is of interest. 

In spars with plywood webs the direction of the fibers is usually 
at 45° to the flanges, since under shearing stress — apart from an 
increase in shear stiffness — a substantial increase in buckling load and 
hence often also in the ultimate load of the web relative to a plate 
fixed parallel can be obtained. 
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In order to include these important practical cases also., the 
problem of buckling of special— orthotropic plates, so exhaustively 
explored by Seydel for shear load, is extended to include the problem 
of buckling of ort hotr op ic plates with arbitrary principal stiffnesses 
relative to the plate edges (termed "general— orthotropic plate" 
hereinafter) . 


II. THEORY OF ELASTICITY OF THE 
GENERAL-ORTHOTROPIC PLATE 


The investigation of buckling of the general— orthotropic plate 
requires a closer insight into the elasticity theory of such plates. 

The general— orthotropic plate is characterized by the inclination of 
the principal stiffness directions (of the principal axes system) of 
the plate with respect to the coordinate system to which the stresses 
and strains are referred. For the special— orthotropic plate, the 
principal system of axes is coincident with the coordinate system 
(figs. 1(a) and 1(b)). The inclination of the principal system of 
axes with respect to the system of coordinates of the plate necessitates 
a construction of the laws of elongation of the general— orthotropic 
plate which, compared to the strain laws of the special-orthotropic 
plate, contains additional terms. Whereas the elastic properties of 
the special— orthotropic plate can be defined by four elasticity 
constants, the Identification of the properties of the general— 
orthotropic plate requires the knowledge of six constants. Hence, 
aside from the stress— strain law — as will be shown in the subsequent 
sections — the differential equation of the bending surface, and the 
equation of energy of the general— orthotropic plate for the bending 
of these plates (both equations being used for the further investi- 
gation), two additional terms are obtained which contain these new 
additive elastic constants. 


1. The General Stress— Strain Law of the Orthotropic 
Plate (Plane Stress Condition) 

(a) The Elongation Equations 

The law of elongation of the general— orthotropic plate can be 
written in the form (reference 5 ) 

e x = ct ii ,CT x + a l2 ,<J y * a l3 , ' r 

e y = a 21* cr x + a 22 t °y + a 23 ' T 

7 = a 31 'a x + a 32 *a y + a 33 »x 


> 


( 1 ) 
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The moduli of elasticity a' for the coordinate system xy are 
dependent on the angle cd of the system* which the principal axes 
system forms with the coordinate system. (These general moduli of 
elasticity are hereafter denoted by a'.^, the special moduli (for the 

case cd = 0) with aq E .) 

It can he proved (reference 5) that the symmetry conditions 
a 12* = ct2p 1 * a i3* = 0031'* and 0-23' = a 32* exist* hence that the 
elongation law of the general— orthotropic plate is represented by six 
elastic constants. 

The moduli a^* and a ' 23 * ^u^ion (l) are a measure for 

the elongations due to shear stresses in the general case (cd f 0) and 
the displacements induced by axial stresses* respectively. 

These moduli a-j^* and “23 1 disappear in the case of -the special 

orthotropic plate (cd = 0) as will be shown later* so that this case is 
characterized by four constants. 

Introduction of the conventional reciprocal expressions for the 
moduli of elasticity 


„ . 1 . „ _ 1 . „ 
11 = E11' ~ = — ’ 


'11 


E 22 ’ -12 “21 Eii 


E 


22 

22 


J a 33 " g 


gives for the case of special orthotropy the equations: 


e x = “ll^x + a 12 a y = ^T( CT x ~ v ll CT y) 


E 11 


e y “ a 12 a x + a 22 a y = ^~( a y “ v 22 a x) 


e 22 


7 = a 33 r = - 


(la) 


On top of that* a-] -| = = — 3 for the case of the isotropic plate* 

E 

so that the equations of elongation read: 


e x = a ll°x + a 12 ff y = | C°x ~ va d 
fy = a i2 cr x + a ll cr y = jjr (o-y — vo^) 


7 = a 33 T = - 


T_ 

G 


(lb) 
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And. since a further relation can he established, between the 
moduli on* a i2* and a 33 > the elastic constants for isotropic 

material are reduced to two. 

(b) Determination of Moduli of Elasticity 

For the general— orthotropic plate, the elastic constants in 
equation (1) are dependent on the angle cu formed by the principal 
axes of the orthotropic plate with the coordinate system. In the 
following, the moduli of elasticity for any system of coordinates 
orientated with respect to the principal axes of the plate are determined 
in relation to the moduli of elasticity for the principal axes direction 
and the angle formed by the coordinate axes relative to the principal 
axes . 


From a general— orthotropic plate element of which one principal 
axis is inclined at angle cd to the x axis, a rectangle is cut out in 
such a way that the edges of the rectangle lie parallel to the principal 
axes of the plate element. The plate element is loaded along its edges 
by the stresses c x , cr and t (fig. 2). These stresses in turn introduce 
stresses crp, a and t d at the cut-out rectangle, which can be 

indicated immediately for the plane stress condition with the aid of 
the conversion formulas 


CT 1 

ct 2 


cr x cos'xd + cry sin'vjo — 2t sin cd cos to 
P P 

cr x sinhn + Oy cos cd + 2t sin cd cos cd 

cr x sin cd cos cd — ay sin cd cos cd + tCcos^cd — sin^oo) 


(2) 


Denoting the diagonal of the rectangular element parallel to the 
x axis with 1, the sides of the rectangle have the lengths sin cd 
and cos 00 (fig. 3(a)). Let Alp and Alg be the displacement 

components of a rectangular point of the element distorted by the 
stresses crp, cs^, and t d along the edges of the rectangle. The 

displacement components of the rectangular point along the x axis, that 
is, the elongation of the rectangular diagonal of magnitude 1, can be 
read from figure 3(L) as 


e x = AL X = Alp cos cd + Alg sin cd 


(3a) 
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The angle of shearing strain through which the connecting line of the 
distorted rectangular point has turned with respect to its original 
position — that is, the x axis — is (fig. 3(13)): 

r - C-ai, sin cd + Alg cos cdJ (3b) 


From a corresponding rectangular element, whose diagonal of length 1 
runs parallel to the y axis, the displacement component of the rectangular 
point and the elongation along the y. axis follow as 


= ALy = AL^ sin cd + Alg cos cd 


(3c) 


The displacement components ALj_ and ALq along the rectangular sides 

consist of displacements due to axial strain, transverse contraction, 
and shearing strain of the rectangular element. 

From figures 4(a) to 4(c), where the several contributions of 
the deformation strain are represented separately, follow immediately: 

(a) The displacement due to longitudinal deformation 


ALf = cos 60 = cos cjd 

Alg = e g sin a) = 02^22 sin cd 

(b) Displacement due to transverse contraction 

ALj = cos CD = tf2 a 12 cos 015 

AL2 = sin cd = cri a i2 s ^- n 60 

(c) Displacement due to shearing deformation 


(4a) 


(4b) 


ai _ y sin sin cd 

A1 i - y 0 — - T o a 33 — ~ 


Al 0 = 7 

2 • o 2 


_ _ ~ COS CD 
- T Q a 33 — g— 


(4c) 


* 

NACA Reviewer’s note: This equation was incorrectly written as 

27 = -Al-j_ sin co + Al 2 cos cd in the original German version of this paper. 
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The expressions (4 a) to (4c) introduced in equation (3a) give 
the elongation <= x 


e x = AL X = Ali 003 03 + ^-2 sln 03 


P P 

€ x = a i a n cos 03 + a 2 a 12 003 03 + T o a 33 


+ 0 ^ 22 . s i n ^> + CT l a l2 + T o a 33 


sin cd cos cd 
2 

sin CD cos cd 


The stresses a-^, and. t q can he replaced hy means of the 
equations (2) hy the stresses cr x , Cy and t , so that the final expression 
for the elongation e x of the general— orthotropic plate reads 

€ x = 0 * [an cos ^d3 + (^ ct 12 + a 33) cos 2 cd + 022 sin^csj 

+ a y^ a ll + a 22) cos2m sin2cjL> + a 12 { cos ^ + ainVi) — 0,33 cos 2 cd sin 2 cDj 
+ t [—2 (a q t cos 2 cd — 022 sin 2 cD) sin cd cos cd 

+ (^2ai2 + a 33) sin 01 cos 03 (^cos 2 cd — sin^cn^J ( 5 a) 

The shearing strain of the general-orthotropic plate follows 
from equations (3h) and (4a) to (4c) as 

7 = o x £—2 (a-n cos 2 cd — 022 sin 2 cn) sin cd cos cd 

+ ^2a^2 + a 33) sin 03 cos cd (cos^cd — sin^cD^J 
+ a y[ _2 ( a H sin^cn — 022 cos 2 cd) sin cd cos cd 
— ^20^2 + a 33) sin 03 cos cd ^cos 2 cd — sin^a^J 

+ t [4 ^a-| + 022 — 2a 12) cos ^03 sin^co + <2,33 ^cos^cd — sin^cs)^ ( 5 b) 


i> 
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The elongation e y can he found by corresponding considerations 
from equation (3c) as , ■ 

6y = q x |(a-| 1 + CL22 ) cos 2 cd sin 2 cD + ( cos^cd + siiAo) — a 33 sin 2 cD cos^coj 

+ 3 y joq_i sin^ao + ( 2 a ^2 + a 33 ) sin 2 CD cos 2 cd + 022 cos^cdJ 
+ t j — 2 (a-i 1 sin^ — 0^2 cos 2 cd) sin cd cos cd 

— (^2a^2 + a 33) cos <*> sin (cos 2 cd — sin 2 cD_)J (5c) 

Comparison of the coefficients of equation (1) and (5) yields for 
the moduli of elasticity of a coordinate system inclined at the angle cd 
relative to one principal axis, the following relations 


a ll* 

= cos^cd + 

( 2a 12 + 

a 33) 

• 2 2 
SimCD COS^CD + 

■ 4 

CX 22 sm td 

CVI 

OJ 

d 

• U 

= 0 tj • ] sirrco + 

( 2a 12 + 

a 33) 

sin 2 cD cos 2 cd + 

h 

OC 22 COSTCO 

H P 

ro 

= (a^ + 

a 22 ) 

cos 2 cd sin 2 cD + 

a 12 G°s 4 a> + 

singes) — a 33 3in 2 cD cos^cd 

a 33* 

- ^(a i;L + a 22 

2012 ) 

cos 2 cd 

sin 2 cD + a 33 ( 

' p 2 ^ 

cos cd — sin CD^Z > 

a 13* 

= ~ 2 (a-11 

0 

COSTCO 

— 022 sin 2 a)^ 

sin cd cos cd 



+ ( 2 cl 12 

+ a 33 

) sin ao 

COS CD 

(cos 2 cd — sin c 

«) 

a 23 ’ 

= -2 (a i;L 

sin 2 cD 

— 022 cos 2 cd) 

sin cd cos cd 



- (2(^2 

+ a 33 

^ COS CD 

sin cd 

(cos 2 cd — sin 2 

cd) 

J 


As essential result of the determination of the moduli of elasticity 
it is found that the six constants of the general— orthotropic plate can 
he derived from the four constants of the special— orthotropic plate 
(principal direction constants an , 0-22 ■> a 12-> 0-33)* 80 "that 

additional experimental determinations of moduli on plates beyond these 
four quantities are not necessary. 
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For isotropic material., the conventional relation between the 
modulus of rigidity. Young’s modulus, and transverse contraction exists 

G = — -= 

2(1 + v) 

It is logical to expect a corresponding relation for orthotropic 
material, so that here also the experimental determination of the 
modulus of rigidity would he superfluous. 

The modulus of rigidity of a general— orthotropic plate is according 
to equation ( 6 ) 

CX33* = + a -22 — ^ a '12) s; '- n ^ a) cos^a) + (X33 ( cos^jj — sin^ai) 2 (6a) 


This reldtion indicates that the modulus of rigidity * is 

independent of the modulus 0,33 valid for the principal direction of 

the orthotropic plate only when u> = •45° (hence, when the directions of 
the principal axes of the plate coincide with the principal stress 
directions in pure shear) and that the dependence is limited to the 
moduli of elasticity and transverse contractions in the principal 
directions. Hence, in this case 


a 33 ' “ a ll + a 22 2a 12 


( 8 ) 


Written in the conventional manner, equation ( 8 ) reads: 


1 _1_ + Vll Vg2 

Eh E 22 e 11 e 22 


E 11 E 22 

E ll(l + Vgg) + e 22 ( e + v ij) 


( 8 a) 


Hertel (reference 7 j p. 135) also derived this equation for this 
particular case (cu = 45°) for the shear modulus of plywood (reference 7 )* 
In all other cases (cu £ 45°) the determination of shear modulus a. 33 ’ 

is predicated on the knowledge of shear modulus 0,33 for the principal 

axis (oj = 0) . Even in the case to = 45° the experimental determination 
of 0,33 cannot he dispensed with, since modulus 0,33 must he known in 

order to determine the moduli of elasticity aq_q’, a 22 1 ’ anc ^ a i 2 ** i- n 
this direction, according to equation ( 6 ). 
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So, in contrast to the isotropic plate, the experimental determi- 
nation of the shear modulus independently of all other principal 

elasticity moduli is necessary in every case. 

In this connection, another particular case of the general— 
orthotropic plate is pointed out: 

If, as may happen in measurements, the relation 


cxii 4* — 2 <x-]_2 


a 33 


exists between the moduli of the principal axes, the moduli of the 
plate in the remaining directions (cd / 0 .) are 




a 12 * 

“ °12 

= Constant 

a 33* 

= a 33 

= Constant 

a ll' 

= a xl 

cosV> + 

( a ll 

<x 22 * 

= a xl 

sin^cD + 

( a ll 

a ll f 

+ 022 ’ 

= a ll + 

022 

a 13 * 

= ag 3 ' 

CJ 

11 

r-f 

d H 

1 


022) sin^cD cos^cd + 022 sin^co 
022) sin^cD cos^od + CC22 cos^cd 
Constant 


Sin CJD COS CD 


( 10 ) 


as is readily apparent from equations (6). 

Thus, in this particular case, the modulus of shear and the 
modulus of transverse contraction are independent of the angle of 
direction, while moduli a/Qt^ 022 * > a,j_ 3 r , an< ^ a 23 ' continue to be 

dependent on cd. 

A particular case of this kind, in which the general— orthotropic 
plate has isotropic properties with respect to a 33 * and a^g’, is 

represented in close approximation by a three-ply beech plywood plate. 
(Compare Section II fig* 18(a).) 

For completeness the derivation of ( 7 ) governing isotropic 
from equations (6) and (10) is included. 


iraterial 



MCA TM 1263 


11 


For isotropic material all moduli must be independent of cd, that 
is, "be constant in every direction. The requirement is met, according 
to equations ( 10 ) when the following relations exist: 


a ll = a 22 

( 11 ) 

“ll + a 22 “ 2a 12 = a 33 

Utilizing the first equation of (11), the second takes the form 

2(<xn — a.^2) = a 33 (12) 


and this equation corresponds, in different form, exactly to the 
relation ( 7 ). 

Equations ( 7 ) show further that the moduli of elasticity 0 ,^ 3 * 
and 0 - 23 * disappear for o> = 0° and 90°, that is^ for the principal 

axes directions. They also disappear for the particular case cd = 4-5°, 
when the principal moduli of elasticity a-| 1 and OC 22 are of equal 
magnitude . 

It is to he noted further that in the case cn = ■4-5° the relations 


a ll' = “ 22 ’ 


and aq_ 3 * = 0 - 23 ' 


are valid. 

(c) The Stress Equations 

To continue the study of the orthotropic plate, the presentation 
of the stress condition in the plate in relation to the strain condition 
is required. These relations can he obtained from equation (l) by 
solution with respect to the stress components. The ensuing equations 
have the form 

a x = a ll* € x + a l 2 * e y + a l 3* 7 

ff y " a 21 * 6 x + a 22 * e y + a 2 3 ^ (13) 

T = a 31 ,e x + a 32 * € y + a 33 * 7 
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The connections between the elasticity moduli a* and the elastic 
parameters a* may be obtained according to the rules of solving 
equation systems. 

Accordingly 


°x = 


A = 


Ai 

— ; tT v 

11 

II 

( 14 ) 

A y 

A 

A 


berminants of 

the system of equations 

(1): 

“11 ’ 

ro 

a 13 * 



o^l* 

a 22 * 

023' 


( 15 ) 

a 3 l' 

a 32* 

a 33 ' 




while the quantities Aj_j and A^ are indicated by the following 
determinants: 



e x a l2' 

a l 3 * 


a ll' 

e x “ 13 ' 


a ll * 

a 12 * 

e 

X 

11 

e y a 22 ’ 

023 » 

A 2 - 

“21 * 

V “ 23 ’ 

a 3 - 

021 ' 

a 22 ' 

e 

y 


7 a 32 « 

a 33* 


a 3i’ 

7 a 33 ' 


a 31 * 

a 32 ' 

7 


(16) 


The solution of the determinants (16), while making use of the 
equations (14), gives 

= ( a 22 ,<1 33' “ a 23*^) e x. + (~ a 12 ,a 33 , + a 13 ,a 23') e y 
+ ^ a i 2 ,a 23 T “ “ 22 * a 13 

^y = (-^12 ,a 33' + a 13 ,a 23') e x + ( a ll ,a 33’ “ a l3 ,2 ) 6 y (IT) 

+ ^a 12 'a 13 » - cx-l! ' 023*^7 

At = ^a 12 «a 23 » - a 22 ,a l3 , ) e x + ( a 12 ,a l 3 ' “ “ll'^'^y 
+ (c^i'022* - a 12 « 2 )7 
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The solution of the determinants gives 


A — *^22 **"*33 * 






•ais’c^s’ 

( 15 a) 


As can "be proved after considerable calculation by introduction 
of (6) in (15a), quantity A — for an orthotropic plate with properties 
defined by the four principal constants — is an invariant, which is 
independent of the position of the coordinate system relative to the 
principal axes, and which gives 


A = 


a- 


*( C 


33\ a ll a 22 “ a 12 


2 ) 


(15b) 


Comparison of the coefficients between the equations (16) and (17) 
gives the connections between the elasticity moduli a* and the elastic 
parameters a': 


^ 11 * 

= ®22 ,a 33' 

[ rr f 2 

- 023 


Aa i2 ' 

= Aa 2 i» = 

-a 12 ,a 33' 

+ ttq 3 ,<x 23 

Aa 22 ’ 

= a 11 »a 33 ' 

' ~ a 13 * 2 


^33' 

= a. 11 *a 22 

' - « 12 t2 


^13* 

= ^ 3 l’ = 

a 12* a 23 ' 

- 

Aa 2 3 * 

= Aa 32 * = 

a 12 ,<x 13 ' 

“ a ll ,a 23 ' 


The equations (18) indicate that the six elastic parameters a* 
of the general-orthotropic plate, like the elasticity moduli a', can 
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be represented with respect to the four elasticity moduli app, a 22 3 
clj_2) and 0,33 of the special— orthotropic plate in the principal 

directions, and which gives, after introduction of equations (6) 
in (18) : 


Ml 


h p p o-li sin^CD + apg cos\d — 2ap2 cos^cd sin^cD 
_ _2_ cos^ai sin^cD + 


33 


a 22* = 


a 12* 


a ll a 22 — a 12 

4 o a ll oos^cd + <X22 singes — 2ap2 cos^cd sin^co 

cos?cd sin^cD + 

p 

a ll a 22 - °12 

(a-pp + ogg^cos^u sin^co — aj^^sirAcD + cos^co) 


33 


4 0 0 

cos cd sin cd + 

33 


a 


33 


a ll a 22 “ a 12 
(sin^tD — cos^cd)^ (ayi + 022 + 2 ap 2 ')co 3 ^tD sin^i 

a 33 


CD 


a ll a 22 ~ a 12 


ap^' = —cos co sin cd 


a 23 ' = cos 03 s ^ ri 03 


2 

<x 33 


/■ 0 0 \ -o-ll sin2cD + a22 cos?cd — ap 2 (sin 2 cD — cos^xd) 

(sin'ncD —cos 4 d ) + — - 

a ll a 22 ~ a 12 2 

r « ~ 

2 / p p \ a-ii cos^cd— cxpp sin^cu— a-ip (sin^ao— Costco) 

(sin cd — cos^cd ) + — — — ^ ^ i 

a, o\ / o 


a ll a 22 ~ a 12 


(18a) 

In the case of the special— or thotr op ic plate, the expressions for a 
are considerably simplified, since ccp^* and Og^* disappear, and result 

in the conventional relations ahready indicated by Huber (reference 5) : 


These two equations were incorrect in the original version of the 
paper and have been corrected by the NACA Reviewer. 
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- _ “22“33 . 

a ll A " 


= 


“22 


E 


’ll 


<x ll a 22 “ a 12 2 1 “ v ll v 22 


“11 


E 


22 


22 A 2 1 _ v V 

“ll“22 “ “12 V 11 V 22 

a - -~ < x 12 a 33 _ ~“12 = V 11 E 22 

12 A a ll a 22 “ “ 12 2 1 “ V 11 V 22 

a ll a 22 “ a 12 2 1 


V 22 E 11 


1 - 


V 11 v 22 


(19) 


*33 


a 


33 


a 13 ” a 23 “ 0 


The atress equations of the special— orthotropic plate read then: 

E, 


J 11 


°y = 


1 ~ v ll v 22 

V 22 E 11 
1 _ v ll v 22 


£x + e . 
1 “ v ll v 22 

E 22 

1 “ V 11 V 22 


€ X + 


■ y 


(20) 


t = G 7 

The further simplified constants for isotropic material read: 


a ll “ a 22 


E 


vE 


2 s 12 ~ 2 
1 - V 1 - 

and the stress equations take the conventional form 


; a 33 - G 


( 21 ) 


Oy = 


E 


1 - V 

E 


( e x + 


a y = 5^2 ( ve * + 


(22) 


T = G7 
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For the subsequent applications, the elastic parameters of a general— 
orthotropic plate for the case cu = 45° are also of interest. 

Equations (l8a) give for this case: 


a n' - a 22 * - “ + r 


1 i a ll + a 22 2a l2 




33 ^ a ll a 22 a 12 


^- = G + t 


11*22 


1 — v v 
- 1 11 22 


a 


12 


i = _ 


1+1 a H + a 22 ~ ga, 12 _ + i E 11 + E 22 + 2V 11 E 22 

a 33 ^ a i 1 a 22 ~ a 12 2 ^ 1 ” 


V 11 v 22 


(23) 


33 


= 1 
4 


“ll + “22 + 2a 12 i E H + E 22 2V 11 E 


22 


l 13 _ a 23 


a ll 0L 22 ~ a 12 
. 1 a ll ~ 


22 


1 V 11 V 22 


! E 22 “ E 11 


a ll a 22 a 12 2 ^ 1 V 11 V 22 


In this particular, case 
a 


11 


a 22 » and 


*13 ' 


= a 


23 


(23a) 


for both the elastic parameters as for the moduli of elasticity. 


It further follows that the constants a-^ ' = a 23* disappear., when 
the principal moduli of elasticity Epj_ and E 2 2 51,8 of equal magnitude. 


The representation of the buckling loads (determined in the second 
part of the report) for the particular case co = 45° and 135 °, respec- 
tively, is referred to the two stiffness ratios. 


D 

D 


23 

22 


t 


t 


a 23* 
a 22 * 


and 



a 12 ' + 2a 33 * 
a 22 r 


A limiting value consideration for these two expressions is added 
at this time. 

This equation was incorrect in the orginal version of this paper 
and has been corrected by the NACA Reviewer. 
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Dpo* 

With equations (23) the stiffness ratios -=^— and . . — i- 5 — can he 

D 22* c 22 f 

represented with respect to the moduli of elasticity of the principal 
axesj they then assume the form 


C 1 - V ^^) + S ( 3 - 2V n) + 3 


^11 


D 22' 


4 C 1 - V 11 V 22 ) + sj C 1 + 2 v ll) + 1 


( 2 *) 


E 


22 


d 23 1 
D 22* 


E 


- 1 


11 


E 


11 


( 3 - 


V 11 V 2 2 ) + 


E 


22 


E 


11 


(l 4 2V U ) + 1 


(24b) 


D 33* 


The limiting value which - can assume for the limiting case E-, 


Doo ■ 

is ^~-r - 3 and for the limiting case Epp = 0 
•°22 


22 


D 


21 


-4 JL_ (1 - v v ^+3 

t E \ 11 22/ J 


D ' 

22 




V 11 V 22) + 1 


The modulus of shear can assume values between 0 and 00 : for 

Dqp* 

G = 0 the result is again = 3 and- for G = °°, =44- = —1- Like— 

D ' Dpp' 


wise, it can be shown that 


22 

D 33* 


22 

D 


33' _ 


= 3 for E 00 = 00 , and = 3 and — 1 

D » dd ’ D 1 

22 ^22 


respectively, for E^p = 0, depending upon whether the modulus G takes 
the value 0 or ®. 
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for 


The- elimination of term 4 -S— ^1 — v-qV. 22 ) from the two equations 

Doq' Dpo' El1 V 

~ 5 ' 3 and — — leaves a relation "between these expressions from 


D 22* 


D 22* 


which the terms with v cancel out: 


D23 * 
4 — - 
D 22* 


•(* 


E, 


22 


D 22/ E 22 


- 1 


E. 


t 1 


11 


The factor 


when 


22 

11 


-‘22 

S 11 


- 1 


E 2 2 


E 


+ 1 1 assumes values "between —1 and 1, 


11 


passes through all possible values between 0 and <*> , so that 


d 23 ’ 

the limiting values of - — - follow the equation 

D 22 


V 

D 1 

22 


1 + 


= ± 


D 22* 


(24c) 


This investigation indicates that the pair of values 
D- - ' 


D 23' 

d 22 * 


and . 33 - exist only in a field that is bounded by the two straight 
D 22 1 

lines 

Dm* 


D 2 3 ' 

d 22 » 


1 + D , 

Ml Md D 23._ 


1 + 


D 33 ’ 


D 22* 


P 22* 


D * 

and by the straight line —33— = 3 (fig. 5a) . 

D * 

22 


4 
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2. Bending of General— Orthotropic Plates 

The theory of "bending of general-orthotropic plates follows the 
classical theory of "bending; hence, first of all, the assumptions that 
the plate thickness is small and that the deflection of the plate itself 
is small compared to plate thickness. 

On the "basis of the coordinate system conformable to figure 5b, the 
intersection moments of the plate read as usual 


ps/2 

M * =J _ a/2 dz; 


ps/2 

— s/2 


V dz; 



TZ dz 


(25) 


The condition of equilibrium between the specific load p of the 
plate applied at an element and the moments M^, My, and M^-y reads: 


* 


2 d2U ™ , ^ 1 

ax 2 ^xby dy 2 


P = 


0 


(26) 


This equilibrium condition is independent of the elastic properties 
of the plate, and holds for the isotropic as for the orthotropic plate. 


The plate is bent by the moments. Assuming, as usual, a linear 
strain distribution over the plate thickness, a point of the plate 
distant by x from the center is, as a result of the curvature of the 
plate, subjected to the distortion 





7 = _2 Z 

dxy 


( 27 ) 


w = deflection of plate . 

* 

Corrected by NACA Reviewer. 
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Introducing the stresses of the general-rorthotropic plate in the 
expressions ( 25 ) for the moments of the plate according to equation ( 15 ) , 
while replacing at the same time the distortions appearing in these 
expressions by the equations ( 27 ) , gives the moments of the plate., when 
the integral is extended over the plate width, as: 



The second derivatives of equations (28), formed and entered in 
the equilibrium condition ( 26 ), give the differential equation of the 
elastic surface of the general-orthotropic plate : 


s3 d^w 

12 r ax 4 


+ 4a 1Q ' - L + ^2a 12 ' + 4a^*) 


13 


dx^ay 


a^w 

a x 2 a y 2 


+ 4a,,- + a t 

3 dx9y3 22 Sy* 


- P = 0 


(29) 
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and, after the abbreviations 


D -* " fi a ll' 


"11 


T) I — 3 q I 

^22 " 12 a 22 


D 


33 ‘ = Ii( a l 2 ' " “33 


+ 2a-,-, 


(29a) 


Ti f — S-3 o f 

13 " 12 13 


= i a 23’ 


the final form of the differential equation as 


11 a x 4 


■w 


a^w .1.^ . a 4 - 

dx 3 ay 


a 4 w 


Dtt ' + 4D 13 f — + 2D_* + to- ' 


33 ax 2 a y 2 


23 


a 4 w 

dxdy~ 


+ D, 


' 2a 57 


= P 


(30) 


Compared to the differential equations of the elastic surface of the 
isotropic and orthotropic plate, respectively, the additive terms with 
the mixed derivatives 

a 4 w „ a 4 w 

and — 

ax 3 a y axa y 3 


appear; they identify the contribution of the elastic parameters a p 3 f 
and ag^’ on the absorption of the load p. 
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For the apecial-orthotropic plate, Dp 3 * and- D 23 1 disappear and 
the remaining constants assume the form indicated hy Huber : 


D-, = Si 


E 


11 


11 12 1 " v ll v 22 


D - s 
v 22 


E 


22 


12 1 - VppV, 


'22 


(31) 


s 3 I E 

2D 33 = ^( v iirT^TT 


22 


E 


+ V, 


11 


22 i __ -y v 

11" 22 V 11 V 22 


+ 4G 


The differential equation of the elastic surface of the special— 
orthotropic plate reads then 


11 i? 


+ 2D 


33 bx 2 dy 2 22 ^ 


V 


= P 


(32) 


For the isotropic plate, when noting that 


Epp — E22 — E 


V - v - V 
11 " 22 


E 


G = 


(1 ~ v) 
2(1 - v 2 ) 


can he put in ( 31 )* it finally gives 


_ s3 E 


Du - D 22 - D 33 - D - i2 i _ 2 
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and the differential equation can be written in the conventional form 


5^w 

57 - 


2d\r 

,2*, ,2 


5\r 


dx <3y^ 


_ P 
D 


(33) 


3- The Internal Flexural Strain Energy of the 
General— Orthotropic Plate 

The internal strain energy of a plate can be generally represented 
by the relation 


2Ai = 


£x e x + Vy + dx dy dz 


(34) 


If the stresses o x , <3j, t in the bracketed term of equation (l) 
are replaced according to equation ( 13 ) by the corresponding distortions, 
the equation for the general— orthotropic plate can also be written in 
the form 


2Ai = 



a ll* e x + a 12 ' e y e x + a 13 ' 7e x 

+ a 2i * 6 x c y + a 22 ' e y 2 + a 23 ' 7e y 
+ a 3l' 7e x + a 32 ,e y 7 + a 33' 7 ^j dx d;7 dz (35) 


If the plate is stressed in bending by moments., the flexural 
distortions e x , e , and y can be expressed again by the corresponding 

« I 

plate curvatures, (equation 21), and gives, after integration with 
respect to plate thickness s, the internal flexural strain energy of the 
general orthotropic plate 



+ 4ag2 * 


a 2 * 

dx 2 


d^w 

dy 2 



d^w d^w 
dx 2 dxy 


(36) 
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For isotropic mater ial, * = a 22* = a H> ^ 4 * & ddition, the 

constants a-^* and a^' disappear and the flexural energy follows as: 



+ 2a 


d^w d^w 



(37) 


If, instead of the constants a^’, a-^, a 33 t ^ Le customary 

reciprocal designations are used again, equation (37) can he replaced 
ty the known expression for the internal flexural strain energy of the 
isotropic plate 



4. Experimental Determination of the Moduli of Elasticity 
of the General— Orthotropic Plate 

The elastic properties of the general-orthotropic plate- are 
characterized by six moduli of elasticity. Inasmuch as, according to 
equation (6), the six moduli of the general— orthotropic plate can be 
determined from the moduli o^, a 22* <x 12-’ an< ^- a 33 f° r ever T direction 

angle 00, the determination of these four moduli is by itself sufficient 
for a complete identification of the elastic properties of the plate. 

But for an experimental check on the relationship between the moduli of 
elasticity represented in equation (6) and the direction angle u> it 
seems appropriate to have a method available which permits the determina- 
tion of all six moduli of the plate. 

(a) Pure Axial Stress Condition — Pure Shearing-Stress Condition 

Nadai (reference 8, pp. 355 ff) pointed out that two stress atti- 
tudes, the case of pure axial bending and the case of pure twisting 
of a plate (figs. 6, 7), are particularly suitable for determining the 
two moduli of elasticity of isotropic plate— shaped bodies. The boundary 
conditions in both cases are easily established and the deflections of 
the plate readily indicated, since the central area of the plate assumes 
a geometrically simple shape in these loading conditions. 
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The determination of the moduli of elasticity from tending and 
torsion tests, especially for plywood plates, which ty reason of their 
inhomogeneous structure over the plate thickness manifest different 
elasticity moduli in "bending than in plane stress, (see section 5 ) has 
the advantage to the extent that the decisive moduli, necessary for the 
investigation of the plate stability, are already obtained this way for 
bending and torsion stress. 

Nadai's method for isotropic plates can also be extended to include 
the general— orthotropic plate and the six moduli of this plate computed 
with it. 

A relation between the stresses introduced by the loading and the 
deflection of the plate is obtained by equating the distortions e x , e 

«7 3 

and 7 according to equations (l) and (27) to each other: 


6 x 2 


7 ~| = “ll'“x + a 12’ CT y + a 13 * T 
dx 


6 = ^ = cl 12 'c x + a^'o + 

J 0 y^ 


(38) -CO 


d 2 w 

7 = -2z = a 13 *a x + 023 + “33 ' T 

dxy 


On assuming a uniform tension stress and shearing stress attitude, 
the stresses in the plate are constant. 


Putting for the stresses in the outer layer z - |j- of the plate 
Oxo > °y = °yo T = T 0 -» e Q. ua ^i o:n (38) can be written as 


Hi ° - I + “l2«y o + “13 '• T o) 

= - f (“12’% + a 22 , °y 0 + “2 3 ' T o) 

d^w If. \ 

~ = - s + “23 CT y 0 + a 33 ’ T o) 


(38a) 
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The general form of a function which satisfies equations ( 38 a), is 



where the remaining terms c^x + Crtf + c^ are disregarded as unessential 

since they merely define the position of the plate in space. Limited to 
the case that the plate is deformed in the plane only by a moment M x , 

it results in a pure longitudinal stress condition and, since a^ Q and t q 
disappear , the equation of the elastic surface — after putting 


xo 


Mg 6 m^ _ 6mx 


W 


bs" 


reads 


w 


_ 6m x 


(cxii’x 2 + a 12 »y 2 + a i3 , xy) 


(40a) 


If the plate is stressed by a moment My only., 


w = (a 12 *x 2 + ^'y 2 + 

s- 5 


a 2 3 ’xy) 


(40b) 


If the plate is stressed by moments only, a XQ = <jy Q = 0, and 

deflection function w is 



a , j c 

23 J 


+ a 


33 


*xy) 


(4oc) 


w 
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On the isotropic plate* a pure longitudinal stress condition* that 
is* a load due to "bending moments or My* produces pure tending, a 

load in pure shear* that is* a load due to twisting moments M^* produces 

a pure twist. (The term with xy in equation (46a) and (40b) disappears 
for the isotropic plate* "because 0 - 13 * = a^' = 0 ** the elastic surface is 

therefore symmetrical to the coordinate axes* that is* pure bending results. 
The terms with x^ and y^ in equation (40c) disappear* the elastic 
surface is antisymmetrical to the coordinate axes* hence pure twist). But 
on the general— orthotropic plate* bending moments* as well as twisting 
moments* even if applied separately* induce bending and twist simultaneously. 

Equations 40a to 40c indicate further that the pure longitudinal 
stress attitude is suitable for determining the moduli 

a-ll'j a 22* 3 “ 12 ' J a 13*-» “ 23 ' 

while the pure shedring stress permits the moduli Oq^'* a 23 *^ a 33 , 
be measured. Consequently* moduli “ 13 ’ and can ^> e obtained by 

either one of the two conditions* but the modulus of rigidity a .33 ’ 

only by the shearing stress, the moduli of elasticity and transverse con- 
traction moduli aqq’* 022 ' ^d “ 12 * on -47 4y means of the pure longi- 

tudinal stress condition. Experimentally* the task can be confined to 
measuring the moduli aqq'* and Oqg 1 for the longitudinal stress* 

and the moduli aq^'* “ 23 ' s-ad a 33 * ^ or ^he shearing stress condition. 

(b) Experimental Determination of the Moduli of Elasticity 

The measurements for the determination of moduli Oqq f , a 22 * * and 
a 12 * a re 3 according to Nadai’s proposal* made on a plate strip supported 

at three points. The bending moment producing the pure lontudinal stress 
in the plate is applied by means of linkages which load the plate in 
three points (fig. 8 ). The deflection is obtained from the difference 
of the readings of three dial gages mounted on the x axis. By 
equation (40a) 


w x = 


_ 6m x 


a 


11 
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for the x axis (y = 0), hence for modulus an when the two outer 
dial gages are spaced at distance x from the central gage. 


W-rS- 5 1 
~ t _ -XL 

a n - 7 0 

0 %; yr 


A corresponding measurement on the y axis gives 


a 


12 


, _ 1 


6m x y 2 


m 


(42) 


while appropriate variation of the direction of the principal axes of 
the plate with respect to the coordinate axes gives 


and 


<x 22 * 


3 

w s J 
1 


6m y 



.“ 12 * 



(43) 


(42) 


respectively. 


The measurements for the determination of the moduli a ^* > 
and a 33 * are made on a square plate which is loaded according to 

figure 9- The loading produces a pure shearing stress since, (refer- 
ence 6, pp. 299 ff) the twisting moments m^y occurring at the 
boundary of the plate can he replaced by two equal and opposite trans- 
verse forces Q normal to the plane of the plate, in the corners of 
the plate. The total single force in a corner is 


P = 2Q = 2m X y ■ 


The measurement of moduli 
a ll * a 12 * * 


a 13* 


023’ 


and 


is similar to that of 
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If the deflection of the plate is measured along the 


(y = 0) 


x axis 


w„ = 


6m xv , 2 
= ^ a 13 x 

3° 


= -3^ (X. *x 2 

s3 


“13 ‘ 


3^c 


3 

2 


(44) 


follows from equation (40c), where x again denotes the distance of the 
two outer gages from the central dial gage. 

Measuring the deflection v along the y axis gives the modulus 

«7 


“ 23 * 


v .V a3 

3Py 2 


(45) 


The modulus of rigidity is defined by the measurement of the 
deflection of the plate along its diagonals (fig. 9). According to 
equation (40c ) , the deflection of the plate in the zero point on the 
diagonal 1 gives — the outer dial gages being mounted at the points 
x = y = b, and x = y = — b 


wp = ( a ]_3 , l :)2 + 023 'b^ + 0,^3 ‘b 2 ^ 

and on diagonal 2 with the corresponding points 1 = b, y = — b, and 
x = -b, y = b: 


Wr 


= 2Z 

a 3 


( a : 


13 


*b + a. 


'23 


— a 


33 


*b 


9 
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Forming the difference leaves 


= 


W 1 “ w 2 ~ a 33 

a 


and hence the modulus of rigidity * as 


. ( W 1 - W 2 ) £ 

a 33 —r " 


6Ph 2 


(46) 


5- Plywood as Orthotropic Plate 

(a) Application of Theory of Elasticity of the Homogeneous General— 
Grthotropic Plate to Plywood 

Since the results of the present report in their practical appli- 
cation are primarily intended for plywood, it is necessary to explore 
the extent to which the theory of elasticity of the homogeneous general— 
orthotropic plate is applicable to plywood. 

Whereas the ordinary wood veneer can, with some Justification, be 
regarded as a plate of orthogonal, anisotropic homogeneous material to 
which the theory of elasticity, and particularly also the theory of 
bending of the orthotropic plate, can be applied, the combination of 
single veneers at 90° to each other into plywood also presents an 
orthotropic plate, but inhomogeneous over the plate thickness. In 
consequence, if the relations applicable to homogeneous material are 
formally retained, the moduli of elasticity of a plywood plate assume, 
in part, different values for stresses in its plane than for stresses 
which tend to deflect the plate out of its plane. 

For a special— orthotropic plate consisting of three identically 
thick - beech veneer plies, for example, whose veneers by themselves show 
the moduli of elasticity E 1 = 180,000 kg/ cm 2 and E 2 = 5,000 kg/ cm 2 

respectively, the assumption of equal strain of the glued laminations 
under longitudinal stresses in the plane of the plate, with transverse 
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contraction disregarded, gives the modulus of elasticity E x in 
x direction (fig. 10) as 




2a. E-, + s. e p 

= sm = 3 1 3 2 


2E-j_ + E2 


F 


= 121,700 kg/cm 2 


and in y direction as 


E y " 


+ 2E 2 


= 63,300 kg/cm^ 


But, stressed in Bending By moments aBout the x or y axis, the 
retention of linear strain distribution over the plate thickness gives 
the moduli of elasticity 



E y « — E-, + -X Ep = 173,500 kg/cm 2 
x 27 x 27 d 

and for E 

v 

E y = ^ El + if E 2 = 11.500 kg/cm 2 


that is, substantially different values than for the plane stress. 
Similarly, the transverse contraction factors in plane stresses are 
different from those for Bending stresses of the plate, as will Be 
shown later. 

The modulus of rigidity of the special— orthotropic plate, on the 
other hand, has, for pure shear in the plane, the same magnitude as for 

* 

These values were altered By NACA Reviewer. 
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the twisting stress of the plate due to twisting moments about the x 
or y axis, since the modulus of rigidity of each lamination is of the 
same magnitude , regardless of whether the grain of the veneer is in 
x direction or at right angles to it (fig. 11). Therefore, the special— 
orthotropic plate is, as far as the modulus of rigidity is concerned, 
homogeneous over the plate thickness. The modulus of rigidity of the 
plywood plate should, in consequence, agree with that of the veneer, but 
on account of the bonding of the laminations the modulus of rigidity of 
the plywood plate compared to the individual veneer, increases according 
to Hertel (reference 7 )j from about G = 7000kg/ cm^ to G = 10000 
to 12000kg/ cm^j depending upon the nature and quality of bonding. 

After the four moduli of elasticity for plane or bending stress of 
the special— or thotr op ic plate have been found — whether by calculation 
from the moduli of elasticity of the veneer (cf . next chapter) or by 
direct strain and bending tests on the plywood plate — the moduli of 
elasticity of the general— orthotropic plate for any angle cd can be 
obtained by the equations (6), so that the theory of the general— 
orthotropic plate remains unrestrictedly applicable also to the plywood 
plate. The extent to which the assumption of linear strain distribution 
is in agreement with actual conditions must be verified by comparison of 
the calculation with the experimental results. 

(b) Effect of Plywood Construction on the Moduli of Elasticity of the 
Special— Grthotropic Plywood Plate 

Investigations concerning the effect of plywood construction on 
the moduli of elasticity of the special— orthotropic plate have been 
made by Hertel (reference 7). Beyond the results of this work, it is 
shown in the following that for the moduli a-| 1 , a 22 , and a^ 2 simple 

approximation formulas can be set up for their dependence on the plate 
construction, and that the formulas derived for the plane stress of the 
plywood plate can also be applied to the plywood plate under bending stress. 

(a) Plane stress of plywood plate .— On decomposing the plywood plate 
into the portion on longitudinally and transversely directed veneers 
(fig. 12), the strain equations for each veneer portion — the moduli,, 
stresses, and distortions of the longitudinally directed plies being 
indicated by superscript (l), those of the transversely directed by’ 
superscript (2) — can be written as follows: 


e x (l) = a ll (l)cr x (l) + a 12 (l) °y (l) 


6 y (l) = a 12 (l)a x (l) + a22 (l) cry (l) 


(^ 7 ) 
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and 


: (2) = a ll (2) °'x (2) + a 12 (2) ffy (2) 


e ^ 2 ) = a i2^ 2 ^ a x^ 2 ^ + a 22^ 2 ^ CT y^ 2 ^ 


and for the plywood plate 


e x “ a ll CT x + “l 2 a y 
6 y = a 12 a x + a 22°y 


( 48 ) 


(49) 


The shearing strain is disregarded, since, according to the fore- 
going, the modulus of rigidity G of the special-orthotropic plate is 
independent of the construction of the plywood. 

From equilibrium, conditions at the plywood plate, with consideration 
of the two veneer portions, follow 


°x ( 1)f 1 + cr x (2) F 2 = a x F 
+ a y ( 2 )F 2 = a y F 


or, when putting. 


E2. = 
FI 


D 


(50) 


•t 1 ) + ex ( 2 )d = a x (l + D) 


* y (l) + a y ( 2 )D = Oy (1 + D) 


( 50 a) 
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On assuming uniform strain over the cross section of the plywood, 
the condition ' 


is applicable. 

Observing that 




a-] 1 (4) 7- ot22^) 


(51) 


oil 


< 2 > . 


22 


and 


a- 


12 


^ , 0,-^1 


12 


equations ( 51 ) in conjunction with ( 50 a), (47), and (48), the stresses 


( 2 ) . ( 2 ) 

a x and v ' 


can he represented with respect to the moduli of elasticity of the veneer 
portion (l) and to the stresses a x and a . The moduli of elasticity 

v 

ct -l 1 and a ^2 follow then from the relation 



a ll a x + a 12 a y - a 22^ 0 x^ + a 12^ 1 ) a y^ 
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The moduli ctgg and likewise a.j _2 can he defined in similar manner 
from the equation 

e y = e y (2) = a 12 a x + a 22 a y = a 12 (l)a x (2) + a ll (l) ^ (2) 


The calculation gives for -| 


a ll 


( 1 ) 


. CL, 

a ll = a H 1 V 1 + D) J 


'22 


Try 


+ D - 


“12 


(l ) 2 


“H^“22^) 


-Dll - 


a ll (l) ' 


“12 


(l) 2 


“22 (1 V 


(1 + D) 


x + Si" 


+ D - 


'“12 


U)> 


a (1) /la (D 

a 22 / \ 22 


la. 


( 1 ) 


(1 + D)‘ 


22 


( 52 ) 


and correspondingly for a^g : 


(1 + D) 


“ll (l) “l2 


(l ) 2 


a 12 = “ 12 ^^ 1 + D ) 


i a, 


„ (U J 


22 a-. 


22 


n “11 (1) |/“ll (l) 

1 + DM + D | 

° 22 (1) A°?2 (1) 


»1 2 (1) ' 2 


1022 


(D> 


(1 + D) 2 


(53) 
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According to Bertel’s measurements on leech veneer, the moduli of 
elasticity a have the following values: 


a ll^ = 5.56 x 10 -6 — 

kg 


022 


( 1 ) = 


200 X 10‘ 


v -6 cm^ 

kg 


= -2.5 X 10 -6 — 
kg 


Since the transverse contraction modulus is sma H compared 

to modulus a 22^^ the terms with MY in equations (52) and (53) 


1 a 


( 1 ) 


22 

can be crossed out, thus leaving the simplified equations 


“11 


( 1 ) 


a- 


“ll 


a (!) 


•12 (1) ( ^ 
+ D " d)„...(i) T - — (i), 


a ll a 22 


a 22 v 


1 + 


dll 


( 1 ) 


• a 22 


( 1 ) 


D 


°11 


( 1 ) 


, a, 


22 


( 1 ) 


+ D 


( 52 a) 


and 


a 12 - a l2 


(^(l + D) 


dll 


a 22 


( 1 ) 

oT 


(1 + D) 


1 + 


an 


a 22 


( 1 ) 

TIT 


a- 


D 


'll 


k a 22 


( 1 ) 

or 


+ D 


(53a) 
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Elimination of the term with 


a i2 (l)2 

a ll (l)a 22 (l) 


In equation ( 52 ) , leaves 


ce- 


ll “ “11 


( 1 ) 


(1 + D) 


1 + 


a ll^ 


a, 


'22 


( 1 ) 


D 


( 52 b) 


The effect which the modulus of transverse contraction exerts on the 
modulus of elasticity of the plywood is quite small and amounts to less 
than 0.5 percent, according to comparisons of the approximation formula ( 52 b) 
with the exact formula (52). The modulus of transverse contraction of 
the plywood is defined by (53 a ) with approximately the same degree of 
accuracy; with the use of (53 a ) and ( 52 b), this modulus can also be 
written as 


V = - 512 = V (!) 1 + 11 

a ll 11 OggC 1 ) 

1 + _££ D 


^l 


( 1 ) 


( 5*0 


The corresponding approximate equations for the moduli ccgg 
a 12 

and Von = — are 

22 0^2 


022 = a 22 


( 1 ) 


1 + D 


CKop(l) 

1 + — D 


( 55 ) 


a li 


(1) 


and 


V 22 - V 22 (1> 


( 56 ) 


1 + 


0t 22 


( 1 ) 
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( g) Bending stress of the plywood plate .— The fact that the tent 
plywood plate haa different moduli from the plane stressed plate can he 
taken into account in simple manner. 


Assuming linear strain distribution over the plate thickness, the 
equations (47) to (49) apply also to the bending strain of the lamina- 
tions, if the elongation of the l ami nations is referred to that with or 
across the grain, (fig. 13 ) (say, to the top boundary layer, for example). 
The elongations of veneer portions in this boundary layer must be equal 
to each other (equation 5l)» The equilibrium conditions for bending 
follow from the consideration that the sum of the portions of the bending 
moment which the veneers support must be equal to the total bending 
moment applied at the plywood plate 


*x (l) W 1 + °x (2>¥ 2 = 

Oy^Wl + O y ( 2 ) W 2 = <Jy¥ 


( 57 ) 


where Wq represents the portions with the grain and W 2 those across 
the grain to the section modulus per unit length (fig. 13 ) 


Vq = 


^s/2 

-s/2 


^ ^(with the grain) 


3 /2 


W 2 = 


's/2 

-s/2 


^ ^"(across the grain) 


V2 


and W is corresponding to the homogeneous material 


W = 31 
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If, as for the plane stress, the factor 



is introduced. 


equation (57) takes the same form as equation ( 50 a) and all the relations 
derived for the plane stress apply to the "bent plate, too; however, it 
should he home in mind that in bending for one and the same plywood 
plate, a different factor D is used for determining the moduli of 
elasticity from that for plane stress. To illustrate: For a plywood 

plate huilt up of three identical veneers in plane stress, D e = 0 . 5 , as 
against Dj-, = O.O385 in bending; for a plywood plate of five identical 
plies, D e = O.667, as against = O.263 (Table 1). 


Figure 14 represents the moduli of elasticity Eqq, Egg, 
transverse contraction factors ^ 22 against the factor D 

which characterizes the construction of the plywood. By this D any 
desired construction of plywood, hence also plywood consisting of plies 
of different thicknesses, can be defined. The D factors for plywood 
of various construction, but identical ply thicknesses (3, 5> 7j 11, 00 plies) 
under plane and flexural stress are shown on the abscissa. The calcula- 
tion for determining the moduli of the plywood was based on Hertel’s test 
data for beech veneers. According to it, the following averages hold 
for the veneer: 


E i;l = 180, 000kg/ cm 2 
E 2 2 = 5 , 000kg/ cm 2 
V 11 = 

V 22 = 0.0125 
G = 7000 kg/ cm 2 


The results of the calculation are Included in Table 1. The comr- 
parison of the moduli obtained by approximation with those obtained by 
the exact equations discloses only minor differences. Table 1 and 
figure l4 also contain the moduli of elasticity for 00 = 45°, the 
determination of which is given elsewhere. 


4o 
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(c) Effect of Angle cu on the Moduli of Elasticity 

Wood veneer and — according to_the arguments of the preceding 
section — plywood as well may be regarded as examples for anisotropic 
(orthotropic) material. Knowing the moduli of elasticity along the 
principal axes, the application of equation ( 6 ) permits the representation 
of the change in the moduli with the variation of angle cd. ' 

(a) Veneer .— Regarding the moduli of elasticity in the principal 
axes directions of beech veneers, Hertel*s test data are available. 

Figure 15 shows the moduli of elasticity a for the veneer with these 
principal axes values plotted against the angle to. 

It is seen that modulus a-^ ' , agg ' > a^g’, a 33 * have a period 

of 90°, while modulus a i 3 * and o^q' repeat only after 180°. The 
values of and a,^* disappear for 0° and 90° (that is, for 

special— orthotropic plates) and agree for co = 45° and 135 ° , respectively. 


The marked anisotropy of beech veneer is particularly apparent on 
moduli a 11 ’ and a^g', which manifest very great -variations with the 

change of u>. 


The moduli an *, cc^g and 0 - 33 1 are shown again in figure 1 6 , but 

in the reciprocal manner of writing Ejj_' = — — — QTlH v - ' = — _ 


“11 ’ 


and V xl t = 


co- 


ll 


The modulus of rigidity decreases within narrow limits and assumes the 


minimum value of G = at a) = 45° • Of interest is the modulus 

cur 

of transverse contraction v-q'j which is positive in the range of the 

principal axes directions, and negative in the range of a) = 11 . 5 ° 
to 78 . 5 °, that is, a veneer strip whose angle cjd lies in this range 
becomes wider under tensile stress, a phenomenon experimentally observed 
by Her tel at 00 = 45° • Besides the measurements of modulus and 

V-j_q at o> = 0°, Hertel also made some measurements of the modulus 

at o> = 45°, which showed E^' = 11,000^^- and V 1 = — O.I 85 . These 

cm 2 

with the computed values of E -p * = 11 , 70 (^ 75 - 

cm^ 


values indicate good agreement 
and v * = - 0 . 17 . 
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(P) Plywood.- With Bertel* a moduli of elasticity for the single 
veneer as basis,, the calculation gives, according to the arguments of 
Section II, 5b, the moduli of the principal axes of the plywood repro- 
duced in table 1, a distinction being necessary between the moduli for 
bending and plane stress. The modulus of rigidity of the plywood is 

increased to 10 , 500 “^ — as a result of the bonding of the laminations, 

cm^ 

as against 7,000^— f or the nonglued veneer. 

cm^ • . 1 


Figures 17 and 18 represent the modulus a* for plywood of various 
construction plotted against angle co, and figures 19 and 20 the corre- 
sponding reciprocal values for plywood plates stressed in plane and in 
bending. 

The equalization of modulus of elasticity Epp and Egg at co = 0 

and 90° as a result of the cross laminations compared to the simple 
veneer is readily apparent; the moduli approach each other with increasing 
number of plies. A balance of the moduli of elasticity in the range 
of 45° angle with the E moduli for 0° and 90° is, however, unobtainable 
on laminations placed at right angles to each other (orthotropic plywood), 
the E modulus drops at cn = 4-5° to about 25 to 30 percent of the 
values Epp at co = 0°, depending upon the number of plies. An equali- 
zation of this range could be obtained also by a plywood which contained 
diagonal laminations in addition to those at right angles to each other. 
But such plywood has found no practical application up to now. 

The modulus of rigidity G for orthotropic plywood shows the typical 

rise from 10,500^2— at cs = 0° and 90° to about 40,000^— at co = 45° 
cm^ cm^ 

The great gain in shearing strength afforded from the use of plywood 
plates having their grain diagonal for shear trusses over normally built- 
up plates is readily apparent. 

The diagrams indicate further that the modulus Epp* for plane 

stress varies within comparatively narrow limits with the change in 
number of plies; in the range of ca = 45° particularly, the variation 
of the Epp* modulus is scarcely noticeable. 

The customary method In practice, independent of the construction 
of the plywood, to figure with a modulus of elasticity of 100,000-^j in 

longitudinal direction (oo = 0 ), of 300 , 000 ^— _ 45 O) p n diagonal 

cm2 


k2 
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direction and of 70 j 00C^S_ in transverse direction (ay = 90°) appears 

cm^ 

to be justified to some extent according to the values of figure 19. 

On the other hand, the moduli of elasticity E - 1 ^ in tending for 

co = 0 and 90° are scattered over a wide range; especially the values 
for three-ply plywood differ considerably from the others. It therefore 
seems no longer permissible to introduce fixed values for the moduli of 
elasticity independent of the construction of the plywood, when plywood 
plates stressed in bending are involved. 

Even the values of G' in plane stressed plywood are not essentially 

dependent on the number of plies, so that here also the value of 40,000^2. 

cm 

used in practice independent of the construction seems Justifiable. 

But the modulus of rigidity in flexurally stressed plywood plates 
is, for a> = ^5° .> very much dependent on the construction of the plate. 

An interesting particular case is represented by the three-ply plywood 
which, independent of the angle ay, has a constant modulus of rigidity 

of G* = 10,500^-. 

cm^ 

The dependence of the moduli G', E-q* = E 2 2 *j v n f = v 22 r for 

diagonally oriented plywood plates (oj = J +5°) on the construction factor D 
is also apparent in figure l4, which shows the minor change of the 
moduli, especially of modulus E’ in plane stressed plates. 


III. THEORY OF BUCKLING OF THE 
GENERAL-ORTHOTROPIC PLATE 

A. EXACT SOLUTION OF THE PROBLEM 
1. The Stability Condition 


If an infinitely long, general— orthotropic, homogeneous and elastic 
plate strip is stressed at its edges by uniformly distributed axial 
forces N x and N y, or shearing forces N^y', respectively, (fig. 21), 
these forces produce components at buckling w of the plate, which are 
at right angles to the plate and of magnitude (reference 6, p. 305): 


-P 


= N„ 


^ + 2Njry — + Ny 

^ dry 7 


d£w 

dy 2 


(59) 
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If this loading of the plate by the external forces (through corre- 
sponding increase of the forces) reaches exactly the value Pj which 
corresponds to the value p as a result of the internal resistance 
against bending according to equation (30), the plate becomes unstable, 
that is, it collapses into corrugations , or wrinkles . 


For this extreme, case of stability, the differential equation of the 
elastic surface becomes, according to equation ( 30 ): 


D-, « £-2- + 4D-, q + 2D„ * — dV .. + * 

dx^ ^ dx^dy dx^dy^ 3 


O^w dK 

— + Dpp’ — j- 

dxdy 5 6y 4 


+ N x + 2N + M & = 0 

X dx 2 7 <*dy y dy 2 


( 60 ) 


The characteristic values of this differential equation give the 
desired critical loads at which the plate begins to become unstable, 
the characteristic functions w corresponding to these characteristic 
values, the buckling forms of the plate. 


2. Formula for Solving the Differential Equation 

For the infinitely long strip, the solution w of the differential 
equation in x direction must be a pure periodic function, in which 
case it is then possible also to apply to the extended differential 
equation of the general— orthotropic plate the Southwell— Skan formula 
(reference 4) for the isotropic plate: 

1*5 iXj (TO) 

w = e a e 


where k is a real quantity. 

■ The formula (TO), introduced in the differential equation, gives the 
characteristic equation 


2iT) * 

l 4 „ u , 3 


Ai + 


+ 


r 2D.no' N a , 

33 K 2__y_h2 + 


d 22 ' 


^22* 


22 


' 4 D 13* 

;d 22 ' 


3 2 V 

c o K 

D 22' 




K _ 


d 22 < 


d 22 » 


= 0 


(Tl) 
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The characteristic equation gives four roots X^, so that for 
given k the complete integral of the differential equation reads: 


w = e 


4 X 

he — 
a 




lx la 


C 2 e 




i\J- 

n 3a , 

C 3 e + 


v 


1X 4 || 


(70a) 


With the roots X^ the characteristic equation can also he written as 
(X - XqMX - X^X - A.3) ( X - X4) = 0 

or resolved 

X^ — X~^(X-^ + Xg + + x,^) + X^~(X-^Xg + X-^X-^ + X-^X^ + XgX^ + XgXij. + X^X]^) 

— X(X.qXgX,^ + x-^XgX^ + X-^X^lj. + XgX^X^) + X-^XgX^X^ = o (72) 

A comparison of the coefficients of the term with X^ of 
equation ( 71 ) and ( 72 ) gives: 


to 23 - 

Xf + Xo + Xo + Xp = — k 

D 22* 

The roots X can, therefore, hfe expressed h j 


(73) 


D 23 f 

X-i o = a ±- 3 — K 

1,2 %2' 


, , V 

X, k = -a ± 7 — k 

3>^ T) | 

u 22 


> 


( 74 ) 
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According to the principles governing the roots of equations with 
real coefficients, a must always he real and ^ and y either real 
or purely imaginary. 

Further comparison of the other coefficients in equations ( 71 ) and 
( 72 ) and introduction of formula (74) gives in place of the character— 
istic equation the relations: 


N a £ 


-2a 2 - 0 2 - y 2 = 2B 1 2 k 2 £ 


D. 


22 


a (p 2 _ y 2) _ 2 B 2 3 k 3 + 


V 2 D 23* + N xy a2 ' 


D 22' ^22* 


J 22 


(75) 


(a 2 -0 2 )(a2- 7 2) = B 4*4 


_ N y »2fe 3 -\ 2 | ^ ya 2 Dg3 . Mxa g 
D 22' \ 22 ) D 22 ' D 22* D 22’ 


The quantities B x , B 2 , and B 3 represent expressions of sums of 
the stiffness quantities of the plate: 


2B-L 2 = 

\ D 22 V 


+ 2 


P 33* 

I> 22 ' 



2 hi 


'22 


f 


I 



( 76 ) 
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3. The Boundary Conditions and Buckling Equations 

As "boundary condition for the plate strip, the two extreme cases of 
the clamped and freely supported edges are investigated. 

In both cases, there is no displacement normal to the , plate at the 
edges. For clamped support, the angle of inclination dw/dy of the 
plate disappears at the edges, while, when freely supported, the fixed 
end moment My must disappear. The boundary conditions for fixed end 
support therefore read 


w = 0 for y = ±a 


— = 0 for y = ±a 


and for free support 


w = 0 for y = ±a 
My = 0 for y = ±a 


(I) 


(II) 


The four boundary conditions (i) and (il) give four linear homogeneous 
equations each, for the determination of the constants C n in equa- 
tion ( 70 a) which, depending upon the boundary conditions, take different 
forms : 

(a) for rigidly restrained support: 


i A.p i 1 2 i A. -n 1 

for w = 0, y = a: Cpe + C 2 e + C^e + C^e = 0 


i A»p i A- 2 iA.q — 1 Aj^ 

for w = 0, y = -a: C^e + C 2 e + C^e + C^e = 0 


1 A .-1 


. '^1 i^2 1^3 

for ^ = 0, y = a: C 1 A. 1 e + CgA^e + C^A-^e + C^A-^e = 0 

w »7 


iA.] 


for $* - 0, , - + c aV - 1Xs + c 3 , 36 - 1X 3 + ^ __ 0 


(KD) 
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Introducing the relation 

e 1 *' n = cos X^ ± i sin X^ 

and equating the real and imaginary parts of (l(l)) to zero results in 
C^cos Xp + C2COS X.2 + C^cos X3 + C^cos X^ = 0 
Cpsin Xp + C2Sin Xg + C^sin X3 + C^sin X^ = 0 


( 1 ( 2 )) 


0 cos X^ CgXgCos Xg "f* C*^X»^cos X^ Gj j Xfj cos X^ — 0 


C-^X-^sin X-^ + CgX^ain 


Xg + C^X^sin X3 + C^X^sin X^ = 0 


These equations have solutions different from zero only when the 
determinant of the coefficients disappears. The determinant put equal 
to zero gives a transcendental equation with the roots X^j the buckling 
equation. 

This buckling equation reads therefore 



cos 

Xi 


cos 

ro 


cos 

X 3 


cos 

x 4 


sin 

Xi 


sin 

Xg 


sin 

X 3 


sin 

x 4 

\± 

cos 

X 1 

^2 

cos 

^2 

X 3 

cos 

x 3 

X4 

cos 


X 1 

sin 

X 1 

X2 

sin 

X2 

x 3 

sin 

x 3 

X4 

sin 

X4 


= 0 (1(3)) 


The solution of the determinant gives 

(Xp — A 3 ) (x »2 — X^)sin(Xp — K^)aln(x^ — X.]^) 

-(Xp - X2)(x 3 - X 4 )sin(Xp - \^)sln (\ 2 - = 0 (l(4) ) 
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Insertion of the expressions (74) for X, n gives, after some 
rearrangements, the buckling equation for fixed end support in the form 

20? (cos 2(3 cos 2 7 — cos 4a) — (4a^ _ p2 _ y2) gin 20 sin 27 = 0 (la) 


(3. Free support 

The first two equations for C n (from w = —0, y = ±a) agree 
with those applicable to fixed end support. 

The third and fourth equations follow from the relation (cf . equa- 
tion 28) . 


a2jj 

“y = D 12* tt + D ^' + 2D 2 3 ’ g = 0 7 


d£w 

dx £ 


d^v 

dr 


(n(D) 


The first term disappears, because 
With the relations from equation (70a) 


= _ &L = 0 for y i a. 
chr a 2 


dfw = _^' l a ^ 
c)y 2 6 


IKf 
a 


4 


n=l 


J n 


4 2 n p iA -l 

■ G 


and 


. x 4 


dxj 


n=l 


n 


v 1 ^ 
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the equations II(l) give for y = ±a the two equations for the 
constants Cq 


SI c n e J®22*^n 2 + 2 D 23 '\n*J = 0 


( 11 ( 2 )) 


— iXjj r~ — I 

>_ C n e h 22 '\n 2 + 2 D 23 «^ K = 0 


Again introducing equation (77) in (ll(2)) and putting the real and 
imaginary parts of the equation equal to zero gives the third and 
fourth equation for On in the form 


c ” oo ° + 2Da3,x "3 ■ 0 

c n sin |® 22 * + 2 D 23 ' A vi^j = 0 


(11(3)) 


Putting for simplification, I>22 ,? 4 i 2 + 2D 23 'A ti k = D n the four linear, 
homogeneous equations for C n in the case of free support read 

Cqcos Aq_ + C 2 C 03 + C 3 COS A 3 + C^cos A^ = 0 

C^sin A-^ + C 2 sin + C^sin A 3 + C^ain A^ = 0 

(II(^)) 

CiDj_cos A^_ + C2l>2 cos ^2 + C3D3COS A3 + C^D^cos A^ = 0 
C 1 D 1 sin A^_ + 02^2 s ^ n ^2 + c 3 D 3 si n A.3 + C^D^sin A^ = 0 


The summation sign was omitted in these two equations in the 
original version of this paper. 


I 
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The determinant of the coefficients of this equation system put 
equal to zero gives again the buckling equation for the case of free 
support. Development of the determinant results in 


(Dp — D3HD2 — Dj + )sin(A,p — A^JsinCA^ ~ I4) 


-(Dp - D 2 )(D 3 - D 4 )sin(lp - \ 3 )sin(l2 - I4) = 0 (ll(5)) 

By (74) the buckling equation for the freely supported plate follows 
after various rearrangements in the form 

{cos 23 cos 27 — cos 4aJ 


— 4 a 2 (p 2 + 7 2 ) — (P 2 — 7 2 ) 2 ] sin 2 p sin 27 = 0 (Ila) 


It is interesting to note that rigidity value D2 3 , / D 22' charac— 

teristical of the general— orthotropic plate cancels out during the 
rearrangements, thus leaving for the general— orthotropic plate the same 
buckling equations from the boundary conditions, as obtained by 
Southwell-Skan for the isotropic plate. 


4. Solution for Pure Axial Compression 

For the case of pure axial compression, the equations (75) take the 
form 


-2a 2 - p 2 - 7 2 


2Bp 2 K 2 


a(p 2 _ 7 2) = 2B ? 3 k3 

(a 2 - p 2 ) (a 2 - 7 2 ) 


- 


n x K 


(75a) 


Nx a 7D 22 


being abbreviated to n 3 
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The three equations (75a) in conjunction with the huclrling' equa- 
tion (la) and (ila) permit the four unknowns, a, p, 7 , and n x , to he 
defined for any, for the present assumed, fixed value of K . The minimum 
value of the values of n^- with respect to « gives the desired 
buckling load of the plate. The stiffness characteristics of the general— 
orthotropic plate enter through the constants Bj_, B 2 , and B^j which, 

in turn, are related to the quantities Dll’, 1 * 22 * ■» D 33 *j ^ 13 * > and -^ 23 * 

in the, qualifying equations ( 75 a). 


From the first two equations of (75 a ), P and 7 are found as 


p 2 = -B-j 2 y £ + B 2 3 ~ - a, 2 

7 2 = — B-, 2 k 2 — Bp^ isl _ a 2 
-l ^ a 


(78) 


and from the third equation follows 


n x = 


B 22 * 


- b 3 V " 


(a- 


p2) (a 2 - 7_) 


(79) 


In the following p is always assumed real, while 7 can become 
real as well as imaginary (reference 1 , p. 280 ) . 

From the transcendental equations (la) and (ila), which for 
imaginary 7 can be written in the form 


2P 


~ Jcos 2 p cosh 2 — — cos 4 aTJ — J^cx 2 — p 2 — 7^J sin 


2 p sinh 2 ^- = 0 (lb) 


and 


8 a 2 p 2 - |oos 2 p cosh 2 ^- — cos 4aj — 

jjla 2 (p 2 + 7 2 ) - (P 2 - 7 2 )^]sin 2 p sinh 2 ^ = 0 


(irb) 


the value a for constant k can be obtained from equation (78) and 
with it the quantity n x from equation (79) . By repeating the 
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calculation for other values of K , n x can he represented as function 
of k , the lowest value of which is the looked— for tuckling load n^. 

For the case of the special— orthotropic plate Fb 2 = 0, = —23., 

1 L d 22 

see equation ( 76)1 a closed solution can he given for the "buckling load 
of the freely supported plate. For the special— orthotropic plate, the 
equations (75 a ) assume the form 


7 2 _ 2B 1 2 fc 2 


2 k 2 
o 


a (p2 _ y 2) = 


(a2 — p2)( a 2 — y 2 ) - — n^tc 2 = -ii — n x K2 


From the second equation of (75^) follows a = 0. The buckling 
equation reads then: 


(p 2 + 7 2 ) 2 sin 2p sinh 2j- = 0 


Since sinh 2 21 > 0, the last equation gives 


sin 2p = 0 and p = g- 


From it, the first equation at (75a) follows as 


,2 . _ 2 2a *2 _ s£ 


« 




The result indicates that the buckling load for pure compression of 
a special— orthotropic plate for ca = 0 and for ca = 90° is of the same 
magnitude, while the buckling lengths for ca = 0 and ca = 90° are of 
different magnitude . 
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The buckling load of the special— orthotropic plate for the case of 
clamped edges^ cannot be given in closed form, but the solution can be 
represented in simple form as function of the "characteristic value" of 


the orthotropic plate 


_ lA ) ll- D 22 
D__ 


and determined in simpler manner 


than in the general— orthotropic case, because a = 0. 


The equations (78) hold for this case also, but the buckling 
equation takes the form 




cos 2p cosh 


1) ' 


a. 


(p 2 + 7 2 ) sin 2p sinh 2 ^- 


Rearrangement results in the transcendental equation 


( 84 ) 


P tan p = - Z. tanh - (85) 

1 i 

from equation (79) follows again 

M x &2 _ Dqi k 2 _ 2 P£ 

D 22 d 22 k 2 

and from (75a) 


2 _ 7 2 4-P 2 

2 231 

D 22 


hence the buckling load 



D 11 D 22 

2D 22 


(7 2 + P 2 ) + 2D 33 


7 2 P 2 

7 2 + p 2 


( 86 ) 


"'"The investigation of the special — orthotropic plate with clamped 
edges does not appear in the literature and is therefore appended for 
the sake of completeness. 
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or 


N x a 2 _ _ d h I) 22 y 2 + P 2 27 g p 2 

D 33 D 33 2 2" 7 2 + p : 


(86a) 


The value N^a 2 /!)^ for every fixed parameter D^D^g/b^ 2 can he 

represented as function of the value 7 regarded as independent variable 
(7 and p are connected through ( 85 ) ) • 

The minimum of this function represents the looked— for buckling 
load. The result of this calculation is Included in table 2 and plotted 
against the "plate value" 


l/ P ll I) 22 = 

D 33 

in figure 22. The closed solution for the case of the supported plate 
is included in the corresponding form in this diagram. 

In Section IIIB the case of the clamped plate which permits no 
closed formal representation of the result is indicated by an approxima- 
tion formula in the form 


which is also included in figure 22 for comparison with the exact values. 
The accuracy of this formula can be substantially enhanced by writing 
(as seen from figure 22) 



N x - 


2n c 

3a £ 


(l/3DilI>22 + °* 88 d 33) 


( 87 a) 


The differences throughout the explored range d = 0.2 to 5 . are 
very small , except in the range d = 0 to 0.2 where they become greater 
and at d = 0 assume a maximum of about 12 percent. At d = 0.2 the 
error is only 12 percent; at d = 5* about 1.8 percent. The values 
between these two characteristic values have an even smaller difference. 
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The case of the Isotropic plate is contained in the special— 
orthotropic plate for the characteristic $ = 1 and gives for the 
supported plate 



and for the clamped plate 


N„ 


17.20D 

a2 


5. Solution for Pure Shear 
In the case of pure shear, the equations (75) read 


-2a 2 - p 2 
a(p 2 ■ 
(a 2 - P 2 ) (a 2 • 
when writing for abbreviation. 


7 2 ) = 
7 2 ) = 

D 22 


2B 1 2 k 2 

2B 2 ^k 3 + n^yic 

b 3 4 k 4 + 2i V A k2 

= n^y and A = 


t 

• 

t 


(75c) 


As in the case of pure compression, the unknown a, p, and 7 can 
be determined from the three equations and the buckling equation (i) 
and (il), when an arbitrary value is assumed for k . After determining 
several k values, the minimum of the function n^y = f ( k) gives the 
desired buckling load of the general— orthotropic plate. From 
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equations ( 75 ) follow after considerable calculations the quantities p 
and 7 


a 2 — 2Aaic + 

a 2 + 2Aa k - jkA 2 K 2 0 L 2 + (4AB 2 3 - + (b x 2 k 2 + 2a^) 2 

(88) 

and, in addition 

nxy = = a (p 2 _ 7 Z) _ 2B 2 3 k 3 (89) 

22 

From the buckling equations (i) and (il), the value a and hence n^y 
is obtained again for chosen values of k with the aid of the equations. 
The minimum value of n^ with respect to K represents the buckling 
load. 

For the case of the Bpecial-orthotropic plate the equations (75b) 
simplify to 


P 2 = ~B]_ 2 ic2 _ 
y 2 = -£ 1 2 k 2 - 


-2a 2 _ p 2 _ 7 2 = 2 ?33 ^ 


22 


a(p 2 - y 2 ) m il^k 


(75d) 


(a 2 - p 2 ) (a 2 - 7 2 ) = -ii. 


22 


because B ^ 2 = ^23., = and A = 0 . 

D 22 5 D 2g3 


The buckling equations are also applicable to this case without 
change. Seydel (reference 3) determined the buckling loads for this 
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case in a very thorough study by means of equations (753.) and plotted 
the result against the plate characteristic. 

For the case of the isotropic plate with the still further simplified 
equations, 

- 2 a 2 - p 2 _ 7 2 _ 2 k 2 

a(3 2 - 7 2 ) = n^K (75e) 

(a 2 - p 2 ) (a 2 - 7 2 ) = 

the solution by Southwell and Skan (reference 4), was obtained with 
Nxv a2 

n^ r = — * — = I 3 .I 0 for the supported plate and n™. = 22.70 for the 
xy p -*-</ 

clamped plate. 


B. THE APPROXIMATE SOLUTION OF THE PROBLEM 


An Increase in the number of elastic constants on the general— 
orthotropic plate compared to the general— orthotropic plate (in pure 
compression the number of elastic constants increases from two to three 
and In pure shear from two to four (see equations (75a) and (75b) 
or ( 75 c) and (753-)) renders the determination of the buckling loads and 
buckling lengths by means of the previously discussed exact method 
extremely time— consuming, since in view of the greater combination 
possibilities of the elastic constants a much greater number of buckling 
load determinations must be carried out. 


1. Method of Solution 

The approximate solution can be obtained by the energy method. On 
transition of the plate from the plane to the buckled equilibrium 
position, hence, on reaching the critical buckling load, the energy 
stored in the plate by the bending strain is exactly equal to the work 
performed by the external forces: 

Ai - A a = 0 (90) 

The internal bending energy owiflfe to the deformation w of the 
plate is indicated by the equations ( 36 ); the energy of the external 
forces during the buckling process may be expressed in the following 
form 
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Equation 10 becomes then 


s3 

12 




013 ^ ^ 


+ 


4^3 • ^ — + ^33 , (rf-) : | 

by 2 bxby \bxby/J 


cbc dy 


-J7 s *( 


'bw'f 

.ox/ 


+ 2 v|| + I v(|) 2 ^ 


(92) 


The condition for stable equilibrium of the plate that the energy 
change during the buckling process shall be a minimum, hence 


BA - b(a ± - A a ) 


can be satisfied by variation of the displacement w of the plate in 
equation (12). This variation process yields the differential equation 
(cf. Ill, A, 1) for the elastic surface of the buckled plate and the 
characteristic values of this differential equation give the exact 
buckling load. But if the buckling load is to be determined only approxi- 
mately, an approximation formula for the buckling form can be introduced 
in equation (12), which differs from the actual form of the buckling 
surface. In this case, a formula is chosen which contains two free values 
which, by application of the condition of minimum energy change during 
the buckling process, can be so determined that the presumed buckling 
form is as close as possible to the actual form. 

A simple formula is applied in the form 


w = A cos sin 2. (x — ey) (93) 

D J-i 

for the freely supported longitudinal edges of the plate, and 

w = ^ ^1 — A cos ^—^sin £ (x - ey) (9*0 

for the fixed end support of the longitudinal edges. 
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This function, originally given "by Timoshenko (reference 6 , p. 365 ) , 
represents infinitely long "buckled surfaces with the half-wave length L 
(fig, 2k) , whose straight nodal lines running diagonally to the longi- 
tudinal axis x form with the y axis the angle cp, the tangent of 
which is equal to ej L and e are the free values. 

Formula (93) satisfies only the geometric "boundary condition w = 0, 
while the "boundary condition My = 0 (see equation Il(l)) is not 

satisfied, "because neither d^w/dy 2 nor d 2 w/cbcdy disappear at the 
"boundaries y = ±a, as a rule. 

On the other hand, formula (9^) satisfies "both w = 0 and — = 0, 

dy 

"but the two formulas still do not satisfy the differential equation of 
the problem. 

The nodal lines of the actual buckling surfaces are, in general, not 
straight lines. But, as even comparatively rough assumptions for the 
buckling form yield satisfactory approximate solutions for the buckling 
load, the simple formulas ( 93 ) and ( 9 ^+) should afford sufficiently exact 
buckling loads. 

As in the exact solution of the differential equation, so the buckling 
load of an infinitely long orthotropic plate strip freely supported, or 
clamped, respectively, will be investigated for pure axial compression 
and pure shear. 


2. Approximate Solution for Pure Axial Compression 
(a) Free Support 

Performing the integrations indicated in the equation (92) by means 
of formula ( 93 ) , the buckling load for the case of pure axial compression 
follows from ( 92 ) as 
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The two free values L 

ditions = 0 and — = rj 
3L 

the conditions — *■ = 0 and 

SL 

SNy 

From — — = 0 follows 


and a are obtained from the minimum con— 
, or, what , amounts to the same thing, from 



dL 


7 


= e 


4 _ k ^ 3 ! 


D, 


22 


•a _ 233' 2 , ^13’ D ll r 

-6^ + 2 — e — 4 — — e + r- 

* D 22 * d oo' Dc 


(96) 


J 22 


from — — = 0 follows 
8 6 



(97) 


Although one unknown L/b can be eliminated from the equations (9 6 ) 
and (97 ) , no closed expression for e with respect to the stiffness 
quantities can be given, on account of the ensuing equation of the sixth 
degree for e; the values e and l/b from equations ( 96 ) and (97) 
must therefore be obtained by a graphical method. 


The introduction of iA/b^ in N x simplifies the expression for 


in 


T>oo’l~ * 3e 2 + - 6 ^2 


* 22 


D 22' 


D 22* 


(98) 


For the special— orthotropic plate, formulas ( 96 ), (97 ) > and (98) 
give the conventional closed relation for the buckling load of an 
orthotropic plate. The terms with D -^ 1 and 1 disappear in 

equation (97 ) > hence e = 0 and equation ( 96 ) therefore gives 


4 
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For N x , after insertion of (99) 

% = 


' 11 D 22 + D 33 


I) 


( 100 ) 


This relation applies exacts since for e = 0, and Dpo’ = Dp., 1 = 0, 

jty jtx ^ 

the formula w = A cos — sin — satisfies the differential equation of 

h L 

the special— orthotropic plate (equation 32). 

A simple relation holds further for the particular case of an 
orthotropic plate of equal principal stiffnesses Dpi* = D 22 *•> 
directions of the fibers of which run at 0° and h5° to the edges, since 


in these cases the terms with 


and ^ 23 ' a 2- so disappear (cf. 


equation (23a)) and 6=0. The stiffnesses D ll' “ 3 22 


and D 33 > 


refer then, of course, to the momentary direction of the axes. 
The buckling load is then 


H x - rf (?22' + D 33') 
b 


( 101 ) 


L _ 


and the half-wave length ^ = 1 


This simple relation does not hold, however, when arranging an 
orthotropic plate of equal principal stiffnesses under arbitrary direc- 
tion of fibers because then e does not disappear. 

Lastly, for the isotropic plate, it gives 


and 


and hence 


D ll* = r) 22 * - D 33* = D 


Di 3 » - d 23 * - 0 


N x = 

D 


■ 


■I II II II ■! 


ii ibii in 11 
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(h ) Clamped Support 


The determination of the corresponding "buckling loads of the plate 
.with, clamped edges in axial compression "by formula (94) gives 




fi 


e + 2D 


33 '(? e 2 + C) _4D23 '(? e3+i,c . 


+ d 22 -i 


e 4 + 8e 2 + 

3 >2, 


( 102 ) 


diL. 

The minimum conditions — ^ = 0 and — — 

bL be 

for defining the free values e and L/h: 


give the two equations 


16 = _ 4 ^ 6 3 + 2 fill e 2 _ 4 V e + ^ 


3 ZS 


D 22* 


d 22 ’ 


D 22 ' D 22 * 


(103) 


'4 t,2 


1,3 3 D 23 ‘ 

2 D 33 * ' 

£ + ^ c 

D 13 * \ 

_ I 

/ v 23 ' \ 

\ D 22* 

D 2 2* 

D 22 7 

1 

( D 22* / 


(104) 


On introducing L^/h^ into the equation for N x j the simplified 
expression 


w x = | ^ D 22’ -U *4 + 3e 2 + - 6 — e 

3 b 


D 22 - 


(105) 


J 22 


is obtained 
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For the case of the special— orthotropic plate, equation (103) gives 


L 1 D H _ 1.316 7 D 11 _ d 65 o 



(10 6 ) 


The buckling load follows from (105) as 


K. 




(107) 


while the expression applies exact for the buckling load of the supported 
orthotropic plate, the expression (107) for the buckling load of the 
clamped plate is only approximately valid, because the formula (97) for 
€ = 0 does not satisfy the differential equation. 

The relation for the buckling load in the particular case of equal 
principal stiffnesses (gd = 0° and ^5°) of the clamped plates reads 


H, 


- ! $ (&11 + D 3 3 ) 


and the half-^wave length 


L = l 

b 2 


jjfr = O.658 


(108) 


The approximate 


expression 


for the isotropic restrained plate is 


8(1 + |/3) *2 8 X 2.731 *2 J2. 

W = UL * D = — % D = 7.28 % D 

X 3 b 2 3 t 2 b 2 


The exact solution is (reference 6, p. 3^5) 
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3 . Approximate Solution for Pure Shear 


(a) Free Support 

As for the pure compression the integrations together with 
formula ( 93 ) give the buckling load 




%L 






(109) 


The expression for the buckling load in shear therefore agrees, up 
to the factor l/ 2 e, with the relation for the buckling load in com- 
pression. Consequently, the equation Sw^/SLe for determining L and e 

must be identical with the corresponding equation ( 96 ). 


The second qualifying equation (iri^/ck gives the relation 


_ L* _ A ^ e 3 + g 6 2 _ ^lll) + 6e 2 _ 2 _ l| = 0 (110) 


2“ 3e -S_ 

L 2 \ • D 22* 


D ' Dxp’ 

22 22 


22 b 


Introducing equation ( 96 ) in (109) again gives the simplified 
expression for N X y 


N. 


Dr 


** _ 7^ 22 


(S- 


D 


3e‘ 


*33 ' . 

- 6 D23 ’ 

) 1 

22 

D 22 V 

Pi 

11 

D 23 = 0 ; hence 


( 111 ) 


N. 


xy 




( 112 ) 
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This equation is identical with Timoshenko’s result (reference 6 , p. 361) 

6 = | ^2 = 0.706 


and 

£ = |/l + a 2 = 1.225 

D f 


hence by (112) 


w = h— «£ -r, _ c C(L -n 

^xy /o 2 D “ 5,66 ~2 D 

l/ 2 b^ b^ 

The exact solution is 

Nxy = 5.35 D, 

the error of the approximate solution amounts to about 6 percent. 

(b) Clamped Support 

The corresponding calculation for the clamped plate gives the 
qualifying equations for e and L; in this instance , equation (103) is 
applicable as in pure axial compression and also the equation 


V-sSal 

■ 2 V D 22 ' 


6 3 + 


D- 


J 22 


e — 


'll 


D 22V 


+ - 2 . 


8 D 33' 


3 D 22 ' 


i6l £ 

3b 2 


= 0 (112) 


The critical shear load is 


N = it Il£_ 
X7 3 eb2 


D 22 4 ^ 3 € 2 + X-^ 

22 V b D 22 * d 22 * 


(113) 
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For the isotropic plate, equation ■ (103) and (112) give e = O. 76 O 
and ^ = 0.826 and hence "by ( 113 ) 


N T = 9.60 


Since the exact solution gives 


N y = 8.98 2 -E 

X o 

b 2 

the error of the approximate solution amounts to 6.7 percent, 


C. BUCKLING LOAD OF THE ORTHOTROPIC ELATE UNDER COMPRESSION AND 
SHEAR LOADING FOR THE PARTICULAR CASE a> = 45° 


In order to minimize the paper work required for the determination 
of the Buckling loads of the general-orthotropic plate, the results for 
several important specific cases are reproduced. 

Of particular interest in practical application is the case of the 
general— orthotropic plate where the directions of the principal axes of 
the plate are inclined at a 45° angle relative to the plate edges. For 
the calculation and especially for the representation of the results, 
this case yields simplifications, since at o> = 45° the relations 
Dj_j_’ = D 22 ’ are applicable (see equation ( 23 a)). The results can then 

Be represented in general form as function of the stiffness ratios D 23 */®!32 * 

and ®33 , / d 22 ,j s i nce J owing to equation ( 23 a), only these two values 

appear in the qualifying equations ((75) and ( 76 )). These quantities 
are determined in the subsequent section (HID) for the special stiffness 
combinations of plywood of varying construction over the entire range of 
angles cu(o> = 0 to l 80 °) in order to demonstrate the characteristic 
relationship between buckling loads, buckling lengths, and angle o>. 


1. Results for Pure Compression 

The buckling loads n^. of the infinitely long general— orthotropic 
plate strip freely supported at the longitudinal edges and stressed in 
pure compression are shown plotted against the parameters Dg^' /^22* 
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D 33 VD 22 * in figure 25 for the specific case oi - 45 ° • The values were 

obtained. by the approximation method described in section IIIB. The 
approximation is made by the graphical determination of the values L/b 
and € from equations (96) and (97) and their subsequent insertion in 
equation (98). 

For purposes of comparison with the results obtained by the exact 
calculation a few of the n^. values obtained by the latter method have 

been included in figure 25. 

The method of determining the exact buckling loads is as follows: 
first, the constants Bp, and B^ for a specific pair of values 

D2 3 , /D22 , j D 33 VD 22 * are determined by equation (76). After these con- 
stants are introduced in equation (76a), the buckling load n x for a 
number of specifically assumed values of K can be obtained with the aid 
of equation (78) and the buckling equation (lib) from equation (79). 

The characteristic behavior of the function ' = f (tc) is shown in 

/ d 23 * Don* 

figure 23 (for the pair — — - = 0 . 5 * 

\ e!2 D 22 

function represents the looked— for buckling load. 

A number of further exact results are obtainable from the qualifying 
equations (75a), without making the previously described calculation, for 
various particular combinations of ^2.7 ^22* ®33 , /®22* : 

(a) When the stiffness ratio ^2.7 1 ^ 2.2 takes the value zero, that 
is, when the plate has the same modulus of elasticity (Epp = E 22 ) along 
both principal axes, the qualifying equations (75 a ) can be written in 
the form 


= v- 


The minimum of this 


.5a 2 — P 2 _ 7 2 . 2 k 2 

d 22 > 


a(|3 2 — 7 2 ) = 0 


( 75 b) 


(a 2 - P 2 )(a 2 - 7 2 ) = - J ~~T ^ ~ n x^ = — n x K 2 


d 22 ‘ 
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The equations agree with the equations ( 75b) for the special 
orthotropic plate, except that instead of the stiffness ratios ^ 33/^22 

and ^h/^ 22 referred to the principal axes directions, the ratios 
■ D 33*/ P 22* D i i */Pp 2 1 appear for o> = 45° • At co = 45 °j D n*/ D 22* 

has the value unity. 


The results obtained for the special-orthotropic plate apply also 


when 


P 23* 
P 22 ’ 


= 0, to the case to = 45° • 


The plate value 


P 33 

(^ 11^22 


hecomes 


then the expression 
figure 25 for 

■ D 22* 


D 33 , / D 22 * • Tlie buckling loads represented in 
= 0 by a straight line are therefore exact values. 


(b) A number of other exact values may be found in simple manner 
for the particular case 



Ml 

p 22 ' 



= 0 


(114) 


The qualifying equations take then the form 
-2a 2 - (3 2 - 7 2 = 2B 1 2 k 2 

a (p2 _ 7 2) = o (115) 

(a 2 - p 2 ) (a 2 - 7 2 ) = 83^^ - 

For the special— orthotropic plate (the buckling equation (la) holds 
in the same form for the general as well as for the special— orthotropic 
plate), the exact relation 


P 22 * 


f (Bq 2 + B 3 2) 


(116) 


is obtained, since 


a = 0. 
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The wave length is 


and 


K = 2. 


JL 

2 Bo 



(117) 


For this particular case between 8X1(1 D 33 , / D 22 , -> the 

relation 


2 




(ll4a) 


Is defined by equation (114) . 
Henue by equations ( 76 ) 



(118) 


Insertion of these values in (116) gives for the buckling load 



and for the buckling length 



(119) 


( 120 ) 
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The equation (119) represents a number of exact buckling loads , 
which in figure 25 lie again on a straight line (denoted with sub— 

D23 * 

script Bp * 0). In the particular case — = 0 and Bo = 0, the 

d 22 J 


values 


n x = * 2 


and 


for 


^33 

B 22' 




valid for the isotropic plate are obtained 


= 0 . 


(c) For the comparison of the exact with the approximate solution, 
it is further of interest that the approximation formula w (equation 93) 
for the buckling surface for the two particular cases treated In (a) 
and (b) satisfy the boundary conditions of the problem and yields the 
buckling loads which agree with the exact values. Accordingly it may 
be presumed that the approximated results in the vicinity of these 
particular cases themselves do not differ very much from the exact values 

Introducing the approximation formula for w 


w = A sin SZ sin — (x — ey) 
,b L 


(93) 


In the boundary conditions for the supported plate 


(1) w = 0 for y = ~ and y = - - 


• + TV. £ 


C^W 




or 


xdy 


( 121 ) 


f or y = k and y = — — 


it follows 

(l) that w = 0 

d^w 


is always satisfied for y = — and y = — 

2 2 


(2) that disappears for y = ^ and y = — ^ leaving 

By 2 2 2 


as 


condition 

D » + 2Dp 9 * & 

22 ^2 23 


* = 0 for y = o 811(1 y ** ~ \ 

ax y 2 2 
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After inserting d^/Sy 2 and d^v/bxdj 3 it is apparent that equa- 
tion (112) can "be satisfied only when the relation 


(123) 


is valid. 


D 2 3* 

The case e = — a2 -f = 0 yields from the approximation equation the 
d 22 


buckling load 


nx . £[1 + ^ 

b V D 2 2 ‘ 


( 101 ) 


already derived in section B, II, 2 and which agrees with the exact 
re suit . 

d 23' 

The relation e = — introduced in the approximate equation ( 96 ) 

D 22’ 

gives the huckling lengths 


Lf. = if, ^33_ 

b 2 2 V D * 

d \ 22 


(120a) 


From equation (97) » further follows the relation 


/® 23 \ ® 33 * 

2 , = -^7 +1 = 0 

\ 22 / D 22* 


and from equation ( 98 ) the buckling load 


(llVb) 


D 33* 

** “ l 3 "w 


(119a) 
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The expressions for L 2 /b 2 and n x obtained from the approximate 

Doo* 

equations with the condition e = ■ a. are therefore in agreement with 

• • - D 22* 

the exact values for the specific case Bp =0. That this actually is 

the case follows from the equation (ll4b) which is identical with the 
qualifying equation Bp = 0. 


These considerations indicate that for Bp = 0 the approximate 

solution agrees with the exact solution as obtained by the differential 
equation. The nodal lines of the buckling surface are actually straight 
lines — which is not the case as a rule — running at an angle with 
respect to the axis, and the tangent of which is exactly equal to 

'23' 


Do-, */Dpp * 


In figure 2 6 , the buckling loads ' n^ are shown plotted against 

^2.^ I ^22* as abscissa and D^’/^p* as parameter. It is seen that 

the buckling load is maximum for — = 0, that is, for plates of equal 

d 22 


stiffness in principal axes directions (cf. equation (24b)), and that 
the buckling load decreases with increasing (absolute) value Dp2 , /^ ) 22 , -» 

that is, increasing difference in stiffness in the principal axes direc- 
tions of the plate, and assumes the value zero in the case of vanishing 

flexural stiffness in one principal axis direction 



The buckling loads lie within an area defined by the range of 
existence of the factors Dp^ ' '/Dpp* and (fig- 5a)* The 

boundary lines of the existence range shown in figures 25 and 2 6 represent 
the boundary lines traced in these diagrams. Since the compressive load 
is symmetrical with the longitudinal!, axis of the plate, the buckling 

loads for oj = 45° ( > Oj and ca = 135° ( < o) are the same. 

\ 22 J \ D 22 * / 


Therefore, the two boundary lines 



1 + 


D 


33 


22 


and 


2gal_ 1 ( 1 + 

D 22 ‘ - - s 



merge into one line. 
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The half-wave length L/b, that la, the spacing of the nodal lines 
of the elastic surface of the buckled plate, was also computed from 
equation (96) by the approximate method and plotted against ^23 '/ D 22 * 

and D33 , /D 2 2 l in figure 27. 

The exact value of the half-wave lengths corresponds to the value k 
obtained in the determination of the buckling loads for which n^. becomei 
a min imu m (fig. 23). The connection between L/b and tc is apparent 
from the equation (70) for the exact formula of the elastic surface; 

j*£. 

since e a is a periodic function, the half wave L must be equal to 
half the period of the function, hence 


L , L 2it 

K — = it, and — = — 
a b k 


(124) 


D 23 ’ 


The half-wave lengths for — — = 0 and B 0 = 0 in figure 27 

n » a 


D 22 

again represent exact values. They indicate that the error introduced 
by the approximate solution is rather small and exceeds no more than 
1 percent in the analyzed cases. 


Figure 27 indicates that for 


D 2 3 t 

D 22' 


0, that is, on plates with 


equal flexural stiffness in both principal axes directions the half- 
wave length is ^ = 1, while with increasing D^' /®22 * * ' tha " t i 3 -> growing 

difference in principal flexural stiffnesses, L/b decreases, and that 
in the extreme case of vanishing flexural stiffness along one principal 


axis 


D23* 

D 22* 


= ±1, L/b drops to zero. The boundary lines of the existence 


range are also shown. 


The value e, that is, the tangent to the angle formed by the 
straight nodal lines of the approximated buckling surface of the plate 
running obliquely to the x axis with the y axis, can be taken from 
figure 28. The straight line according to the approximation equation (93) 
reproduces in general, the actually forming nodal line only approximately. 

D 23* 

It is only in the cases = 0 and B^ = 0 that the nodal lines at 


D 22* 


angle e = 0 and e = 


D 23 ? 

d 22 ’ 


are actually straight lines, because in these 


instances the approximation formula for the elastic surface is exactly correc 
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Figure 28 indicates that the value e becomes zero also for 
= 0 , that is, the nodal lines are at right angles to the edges. 


as 


is readily apparent in view of the symmetrical structure of the plate. 

In all other cases, the nodal lines run obliquely to the edges even 
under pure compression, similar to the case of the buckled isotropic 
plate loaded in shear. As •^23 , /®22* increases, the value e increases. 


Do <3 ■ 

to reach in the extreme case of ■ = ±1 the value unity, the slope of 

® 22 ' 

the nodal lines to reach the value 45°. The boundary lines of the 
existence range are included. 


2. Results for Pure Shear 

The buckling loads 11 ^ of the infinitely long orthotropic plate 

strip, freely supported at the longitudinal edges and stressed in pure 
shear, are represented in figures 29 and 30 for to = 1+5° and co = 135 ° 
plotted against the stiffness ratios ^> 23 ' /^ 22 .' and D 33 , / D 22'* T ^ e 

results were obtained by the approximate method. 

To estimate the error of the approximate relative to the exact 
calculation, the buckling loads for a number of stiffness combinations 
were computed by the exact method and plotted in figure 29a. In table 3* 
the exact result Is compared with the approximate values. The com- 
parison indicates that the result of the approximate calculation, espe- 
cially at higher values of is no longer in as good an agree- 

ment with the exact values as for compressive loading. In the particular 
d 23 ' 

case — — = 0, it is found that the exact results for the special— 
d 22 

orthotropic plate can equally be applied to the plate oriented at 45°« 

For the case of pure shear with / D 2 . 2 ' > the qualifying equa- 

tions (75c) take the form 

-2a 2 - p 2 - 7 2 

a (p2 _ 7 2) 

(a 2 - P 2 )(a 2 - 7 2 ) 


D 33* 2 

D 22* 


= %y K 
= 
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The equations are in agreement with equation (75d) of the special— 
orthotropic plate , when hearing in mind that, instead of the stiffness 
ratios ^ 33/^22 ^ 11^22 re f erre(i to the principal axes directions, 

the values and D-| '/Dpp* referred to cd = 45° appear now, 

D ni * 

and that = 1 for 00 = 45°. 

d 22 * 


The exact buckling loads computed for the special— orthotropic plate 

Doo * 

by Seydel (reference 3 ) are also valid, provided — -2— = 0, for the 

D 22* 

orthotropic plate with direction of fibers at 45°. The plate value 
D 33 //D 11 D 22 a g airi changes to D 33 * /- D 22 l • 


The unusual feature in figure 30 is the fact that the magnitude of 
the buckling loads in shear is essentially dependent upon the position 
of the principal axis of the greater bending stiffness with respect to 
the direction of the shear loading. Figure 31 represents two orthotropic 
plates of equal stiffness, one with the principal axis of maximum 
stiffness at 45°, the other at 135 ° relative to the plate edges. 


D 23* 

At a) = 45°, E 11 > e 22 , as seen in figure 31, and — - f - < 0 according 


22 

to equation (24b); the buckling load — for equal stiffness ratio 


p 33 f 

P 22* 


assumes, according to figure 30 , substantially higher values than for 


CD = 135 °, where < E 22 and hence 


<23 

D 22’ 


> 0. 


In practical application, it means that by suitable orientation of 
the principal axis of maximum stiffness relative to the direction of 
the shear loads a substantial increase in the critical shear load of 
the plate can be obtained. 

Figure 31 shows the ratio of the shear loads for these two cases 
plotted against the stiffness ratio 1 ^* 22 * ' T ^ e i ncrease i n buckling 

load of one over the other is readily seen. It assumes very considerable 

D2 3 « 

values with increasing — that is, with increasing difference of 

D 22 

stiffness in the directions of the principal axes. 
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In figure 32 the approximated half-wave lengths of the "buckling 
surface formed under the critical load are represented in comparison 
with the exact values for several stiffness combinations of table 3 - 
The agreement is very good. 

The approximate value e of the slope of the nodal lines of the 
"buckling sUjTface can "be taken from figure 33* 


3. Buckling Loads of Plywood 

The stiffness values for o> = 45° were computed on several plywood 
plates of different construction with the moduli of elasticity indicated 
In figure 14, and the buckling stresses, half-wave lengths and nodal 
line inclinations determined for compression and shear loading with the 
aid of figures 25 to 28 and 29 to 33* The results were then plotted 
against the plywood construction In figures 3^ and 35* 


4. Example 

The application of the results described in sections 1 and 2 is 
illustrated on a worked— out model problem: 

Consider a very long, freely supported plywood plate built up of 
three identically thick plies, the grains of which are oriented at 
03 = 45° and 135° relative to the plate edge . The plate is b = 2a = 10 cm 
in height, by s = 1 mm in wall thickness. 

Young's modulus E for the principal axes of the plate can be 
determined by experiment or taken from figure 14. 

According to figure 14, they are: 


Eh = 176,500 =2- 

cmr 

e 22 = 11,500 

cmr 


V-Q = 0.196 

V 22 = 0.0130 

G = 10,500 — 
cm 2 
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With these data, the stiffness values for the axes at 45° can he 
computed by (29a) and (33) 


t. _ s- 


Dll f - D22 * = — l G + — 


1 E 11 + E 22 + 2v 11 E 22 


4 1 - 


V 11 V 22 


3 *5 

~ (10,500 + 47,300) = g 57,800 = 4.81 kgcm 


J 33 


-G + p 


1 E 11 + E 22 + 2v 11 E 22 -l E 11 + E 22 “ 2V 11 E 22 


1 — V V 

11 22 


+ — ' ' 

1 “ V 11 V 22 


3 3 

S / n r\ ar\r\ ■ 1. r-r r\r\ . aaa \ S ^ 


= (-10,500 + 47,300 + 91,000) = £- 127,800 = 10.64 kgcm 

12 12 


I) 13 , = d 23' = }§ ( 4 — — ) = “ f| 4 °.550 = -3.38 kgcm for o> 


^ “ v ll v 22 


= Sf. 40,550 = 3.38 kgcm for oo = 135° 
12 

With these values the stiffness ratios become 


5H3_ = ± 3^38 = ±0>T oi 
D 22 4.81 


l33_ = 10^ = 2.20 

D ' 4.81 

22 
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The buckling load of the plate stressed in compression for this 
pair of values follows then from figure 25 as 


_ 

■“x r, 1 
^22 


= 6.89 


and hence the buckling stress as 


_ N x _ n x D 2 2' _ 6.89 X 4.8l 


CT x = -f = - V = = !3.3 tg/cm 2 

a^s 25 x 0.1 


The half— wave length of the wrinkles follows from figure 2J: 


L = - b = 0.817 x 10 = 8.17 cm 
b 


For the plate stressed in shear., figures 29 and 32 show 


(1) for co = 45° 


r = = ^.y^gg ' = 4.96 x 4.81 = 


kr 


3 a 2 s 25 x 0.1 


8.8 kg/cm^ 


and 


L = b = I.76 X 10 = 17.6 cm 
b 


(2) for cd = 135 C 


t,_. = V^ L* = 31.-65. X .1.-81 , 61 . 0 Wcm 2 


kr p 

a^s 


25 X 0.1 


and 


L = b = 0.89 X 10 = 8.9 cm 

D 
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D. BUCKLING LOADS OF THE GENERAL ORTHOTROPIC PLATE REFEREED TO 00 
( PLYWOOD OF VARYING CONSTRUCTION) 


In view of the large amount of paper work involved in a general 
representation of the buckling loads of the general orthotropic plate 
relative to the stiffness factors and angle co, due to the great number 
of parameters j the characteristic relationship between the buckling 
loads and the angle co is exemplified on a worked out model problem. 

The results were obtained by approximation method., the known exact 
values for co = 0°, 45° , and 90° are also given for comparison. The 
results are again limited to a very long plate strip freely supported 
at the edges. 

The calculation is based on the modulus of elasticity of plywood 
stressed in bending, indicated in table 1. 


1. Results for Pure Compression 

In figure 36 the buckling stresses in pure compression are repre- 
sented for co = 0° to 180°. Owing to the symmetry of the compressive 
loading, the buckling stresses are symmetrical to the value co = 90°. 

They are of equal magnitude for plywood with and across the grain 
(co = 0° and 90°) and increase for intermediate angles co. The maximum 
buckling stress for plywood of very many laminations (00) is reached at 
co = 45° and 135 °j respectively. 

The half-wave lengths of the plate represented in figure 37 indicate 
that, if the plate consists of a few plies, the buckling length, with 
the grain, is greater than across the grain. The buckling load increased 
by the greater longitudinal stiffness at co = 0° relative to co = 90°, 
is reduced again by the greater buckling length so that the buckling 
loads are the same again in spite of the substantially different longi- 
tudinal stiffnesses for co = 0° and co = 90°. 

The ensuing average nodal lines which, as under compressive loading, 
are inclined to the plate edges according to section IIIc, and at right 
angles to the edges, in a few specific cases only, are shown in figure 38* 


2. Results for Pure Shear Loading 

On account of the unsymmetry of the shear loading the buckling 
stresses are no longer symmetrical to co = 90°. 



NACA TM 1263 


81 


The close relationship existing between shearing stresses and 
angle to — especially for plywood with few laminations — is seen from 
figure 39. Owing to the unsyrametry of the plywood structure as well as 
to the shear loading, the maximum - shearing stress occurs at about 60 ° 
rather than for od = 45°. 

The half-wave lengths and the slope of the nodal lines in fig- 
ures 40 and 4l also exhibit a marked dependence on angle the half- 

wave length — for three— plywood — fluctuates between ^ = 2.64 and O. 56 . 

For an estimate of the error introduced by the approximate as 
against the exact solution, the exact values for to = 0 ° and 90 ° are 
represented in figures 39 and 40 according to the results by Seydel and 
for ca = 45° according to figure 29. For cu = 45° the maximum error 
is 12 percent; for other cases the error is less. The agreement between 
the approximated and the exact half-wave lengths is remarkably good. The 
approximate calculation which involves incomparably less paper work than 
the exact solution yields a very good insight into the relationship 
between the buckling stresses and buckling lengths and angle o>, notwith- 
standing the cited differences. 


Translated by J. Vanier 
National Advisory Committee 
for Aeronautics 
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Plate Stressed in plane 


Plate stressed in lending 


Number of 
plies 

D 

0 

0 

11 

3 

co = 45° 

D 

' 

O 

O 

II 

3 

a> = 45° 

% 

E 22 

V 11 

V 22 

G 

*11 ' = E 22’ 

■ 

O' 

Ell 

e 22 

Vll 

v 22 

G 

E ll' = E 22' 


G' 

B 

■ 

180,000 

5,000 

0.45 

0.0125 

10,500 

11,700 

-0.113 

4,750 

0 

180,000 

5,000 

0.45 

0.0125 

10,500 

11,700 

-0.113 

^,750 

3 

a 

121,600 

63,400 

0.0355 

0.0185 

10,500 

33,780 

0.604 

40,670 

0.0385 

O 

O 

lf\ 

m 

t— 

1 — 1 

11,500 

0.1959 

0.0130 

10,500 

21,570 

0.035 

10,520 

5 

0.667 

110,000 



75,000 

0.0301 

0.0205 

10,500 

34,210 

0.626 

43,520 

0.263 

143,600 

41,400 

0.0544 

0.0157 

10,500 

31,900 

0.516 

31,370 

7 

0.75 

M 

O 

VJl 

O 

O 

O 

80,000 

0.0281 

0.0214 

10,500 

3^,330 

O.63I 

44,340 

0.4o6 

129,500 

55,500 

0.0405 

0.0174 

10,500 

33,300 

0.576 

37,940 

U 

0.833 

100,400 

84,600 

0.0266 

0.0224 

10,500 

34,400 

O.635 

44,820 

0.573 

116,200 

68,800 

0.0327 

0.0194 

10,500 

34,050 

O.616 

42,180 

CO 

1.00 

92,500 

92,500 

0.0243 

0.0243 

10,500 

34,460 

0.637 

45,140 

1.000 

92,500 

92,500 

0.0243 

0.0243 

10,500 

34,460 

0.637 

45,140 
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TABLE 2.- BUCKLING LOADS OF THE SPECIAL— ORTHOTROPIC PLATE STRIP, 
IN PURE COMPRESSION VERSUS CHARACTERISTIC VALUE OF THE PLATE 


Characteristic 
value i9 


Supported 

plate 


Clamped plate 

approximate ion according to 



0.1 


0.2 


0.3 


0.4 


0.6 


0.707 


0.8 


0.9 


1.0 


1.414 


1.582 


2.0 


2.5 


3-33 


5.0 


Exact values 

Exact values 

Equation (87) 

~ = 4.935 

y- = 4.935 

6.6 0 

5.43 

6.69 

7.74 

5.92 

7-98 

8.88 

6.42 

9.15 

10.02 

6.91 

10.35 

11.16 

7.90 

12.64 

13.44 

8.42 

13.88 

14.58 

8.88 

14.92 

15.72 

9.38 

16.06 

16.86 

:t 2 = 9.87 

17.20 

I8.00 

11.91 

21.96 

22.72 

12.74 

23.95 

24.61 

14.81 

28.52 

•29.40 

17.27 

34.00 

35.10 

21.38 

43.22 

44.58 

29.61 

61.75 

63.60 


5.80 


6.94 


8.08 


9.22 


10.36 


13.88 


14.92 


16.06 


17.20 


21.92 


23.82 


28.60 


34.30 


43.80 


62.80 













































































TABLE 3 .- COMPARISON OF THE RESULTS FOR PURE SHEAR (<d = 45 °) ACCORDING TO THE 
APPROXIMATE METHOD WITH THE ACTUAL RESULTS 


1 

B 1 

^ 3 _ = l.O 

D 22 ' 

v 

— = 0.5 
D 22’ 

B 1 

V 23 _ 

D 22 ' 

— 

V 

= - 0.5 

D 22' 

V 

= - 1.0 

D 22' 

d 33 * 

%>' 

Exact 

Approximate 

Exact 

Approximate 

-*-Exact 

. 

Approximate 

Exact 

Approximate 

Exact 

Approximate 


L 

b 

n XJ 

L 

b 


L. 

b 

V 

L 

b 


L 

b 

”xy 

L 

b 

"xy 

L 

b 

V 

L 

b 

V 

L 

b 

“xy 

: i 



L 

b 

0 

I 

B 

B 

Bfl 


gj 

Qj 

B 

8.125 

1.025 

8.77 

1.021 

g 

■ 

■ 

— 

gM 

■ 

mm 

■ 


B 

1 


1 I 

| 


B 

n 

I 

B 

a 

Efl 

2.32 

M 



13.165 

1.245 

13.98 

1.225 

20.96 

1.45 

23.40 

1.45 


Ml 


I I 

B B 

m 

B 

1 1 

2 

B 

fl 

a 

B 

8.07 

1.126 



17.25 

1.462 

18.22 

1.455 

24.87 

1.69 

27.40 

1.66 

g 

■ 

■ 

■ ■ 





3 

0 

0 

0 

0 

12.06 

1.45 

12.35 

1.414 

20.50 

1.715 

21.80 

1.682 

28.19 

B 

30.78 

1.86 

34.56 

1.98 

39.50 

2.00 


'''According to Seydel (reference 3), figure 4 a and 6a. 
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Figure 5(a).- Range of existence of factors 
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Figure 8.- 


Experimental setup for defining modulus of elasticity 

and a 12 '. 


a ll' » a 22' » 



Figure 9.- Experimental setup for defining modulus of elasticity a-^g', cogg' , 

and cogg®. 









Figure 13.- Assumed flexural strain distribution of plywood 
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Figure 15.- Relationship between modulus of elasticity a' of a veneer and w. 







Figure 18(a). 



Figure 18(b). 
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Relationship between modulus E-q' and G' 
3 plywood plates and w . 
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Figure 19.- 


of plane stressed 


NACA TM 1263 


101 



Figure 20.- Relationship between modulus E-^' and G' of flexurally 

stressed plywood plates and w. 
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Figure 24.- Approximate assumptions for the buckling surface of the general-orthotropic plate for freely 

supported and clamped plate edges. 
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Figure 25.- Buckling load of orthotropic plate strip under axial forces 
(u = 45° and 135°) (approximate values for free support). 


I 







Figure 26.- Buckling load of an orthotropic plate strip under axial forces 
(w = 45° and 135°) (approximate values for free support). 


L 
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Figure 27.- Half-wave lengths of an orthotropic plate strip under axial forces 
(w = 45° and 135°) (approximate values for free support). 
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Figure 28.- 
(u 


Inclination of nodal lines of an orthotropic plate under axial force 
= 45° and 135°) (approximate values for free support). 
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Figure 29(a).- Buckling load of orthotropic plate strip under shear load 
( co = 45 0 and 135°) - exact values for free support. 



Figure 30.- Buckling load of orthotropic plate strip under shearing force ( u = 45° and 135°) (approximate 

values for free support) . 
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Figure 31.- Relationship between the buckling load of an orthotropic plate strip 
(w = 45° and 135°) and the position of the direction of the maximum 
flexural stiffness to the direction of the shearing forces (approximate values 
for free support). 




Figure 32.- Half-wave lengths of an orthotropic plate strip loaded in shear (w = 45° and 135°) (approximate 

values for free support). 
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Figure 33.- Inclination of nodal lines of an 
orthotropic plate strip loaded in shear 
(u = 45° and 135°) (approximate values 
for free support) . 
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Figure 34.- Buckling stress, buckling length 
and nodal line slope of a freely supported 
plywood strip loaded in compression 
(u = 45° and 135°) plotted against the 
construction of the plywood (approximate 
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Figure 35.- Buckling stress, buckling length, and nodal line slope of a freely- 
supported plywood strip loaded in shear ( w = 45° and 135°) plotted 
against the construction of the plywood. 
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Figure 36.- Buckling stress of freely supported plywood strip loaded in 
compression plotted against construction of plywood and u (approximate 
values ) . 
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Figure 37.- Buckling length of freely supported plywood strip loaded in 
compression plotted against construction of plywood and w (approximate 
values) . 





Figure 39.- Buckling stress of freely supported plywood strip under, shear 
load plotted against construction of plywood and w° (approximate values). 
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Figure 40.- Buckling length of freely supported plywood strip loaded in 

shear plotted against construction of plywood and w° (approximate values). 
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